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ABSTRACT 


Synoptic  Interplanetary  scintillation  (IPS)  observations  were 
taken  during  the  summer  of  1976  and  authumn  of  1977  on  the  University 
of  Iowa  COCOA-Cross  radio  telescope  (34.3  MHz),  with  supplementary 
observations  from  the  University  of  Maryland  TPT  array  (38  Miz).  A new 
high  sampling  rate  (10  s *)  digital  system  made  It  possible  to  recon- 
struct the  IPS  power  spectrum  between  0. 1-3.0  Hz.  The  observations, 
combined  with  earlier  (1974)  measurements  of  Integrated  IPS  power 
(scintillation  Index),  have  lead  to  the  conclusion  (based  on  theo- 
retical modelling)  that  prediction  of  geomagnetic  activity  and  associ- 
ated variations  in  energetic  solar  particle  events  Is  feasible  with  a 
lead  time  of  about  24  hours.  The  technique  depends  on  the  observed 
broadening  of  the  IPS  power  spectrum  as  solar  wind  density  enhancements 
approach  the  earth.  This  effect  has  been  documented  for  both  co-rotating 
and  solar  flare-associated  plasma  disturbances.  Full  validation  of  the 
prediction  technique  was  not  possible  due  to  extensive  down-time  caused 
by  unusually  adverse  weather  conditions  at  the  Clark  Lake  Radio  Observa- 
tory In  both  1976  and  1977.  Consequently  about  half  the  available 
resources,  fiscal  and  manpower,  had  to  be  devoted  to  a major  upgrading 
of  the  COCOA-Cross  array.  In  its  improved  condition,  and  using  new 
techniques  developed  from  the  observational  and  theoretical  work  per- 
formed in  1976-1977,  the  COCOA-Cross  may  be  capable  of  extending  the 
prediction  lead  time  for  geomagnetic  disturbances  significantly  beyond 
24  hours. 
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We  proposed  a 24-month  Investigation  of  the  feasibility  of  using 
interplanetary  Intensity  scintillation  (IPS)  observations  of  natural  celestial 
radio-sources  at  34.3  MHz  for  the  prediction  of  solar  particle  events  and 
geomagnetic  activity.  The  IPS  technique  is  the  only  means  available  for 
locating,  mapping  and  tracking  large  scale  disturbances  In  the  solar  wind 
which  give  rise  to  significant  geophysical  activity;  thus  It  can  potentially 
serve  to  provide  continuity  in  the  hierarchy  of  solar-terrestrial  observations 
which  began  at  and  in  the  immediate  vicinity  of  the  sun  itself  with  a variety 
of  radio,  optical,  ultraviolet  and  x-ray  measurements,  and  end  at  or  near  the 
earth  with  near-earth  space  probe  measurements  and  various  types  of  geophysical 
data. 

The  COCQA-Cross  Radio  Telescope,  with  which  we  proposed  to  carry 
out  the  IPS  observations,  was  designed  specifically  to  exploit  the  IPS 
technique,  i.e.  the  collecting  area  is  large  enough  to  allow  measurement  of 
IPS  parameters  on  a large  number  of  radio  sources,  the  beam  pointing  is 
electronic  to  allow  re-pointing  to  a large  number  of  sources  each  day  without 
spending  an  inordinate  amount  of  time  slewing  (which  would  be  required  for  a 
dish  type  instrument),  and  post-detection  electronics  are  designed  to  measure 
IPS  parameters. 

At  the  time  that  we  proposed  to  undertake  the  feasibility  study, 
the  COCOA-Cross  had  demonstrated  both  its  operational  and  scientific  capa- 
bilities, i.e.  more  than  80  sources  had  displayed  IPS  activity  and  daily 
observations  extending  over  6 months  had  revealed  a succession  of  13  recur- 
rent IPS  activity  enhancements  with  a recurrence  period  of  27  days  ( Erekine , 
et  at.,  1976)  which  were  well  correlated  with  geomagnetic  activity. 

Other  work  in  mathematical  modeling  of  radio  scintillation 

(Mitchell  and  Roelof,  1976)  had  demonstrated  that  power  spectra  of  the 

intensity  fluctuations  were  necessary  for  the  proper  interpretation  and 

complete  understanding  of  the  intensity  fluctuation  data.  Until  this 

study  began,  the  data  collection  capabilities  of  the  COCOA-Cross  were 

restricted  to  recording  the  total  integrated  Intensity  fluctuation  power 
2 

m , where  m is  the  scintillation  index.  Therefore,  in  order  to  conduct 
the  feasibility  study,  the  group  at  JHU/APL  (principally  B.  L.  Go two Is) 
constructed  a data  acquisition  system  capable  of  recording  high  time 
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resolution  data  on  magnetic  tape  both  from  the  COCOA-Cross  and,  when 
observing  time  was  available,  simultaneously  from  both  channels  of  the 
University  of  Maryland  TFT  radio  telescope  (also  at  Clark  Lake)  at  38  MHz. 

The  capability  of  collecting  data  on  the  TPT  was  originally 
included  as  a cross-check  of  the  two  instruments,  and  as  an  Improvement 
of  the  statistical  stability  of  the  observations  on  the  limited  number 
of  sources  which  the  TPT  was  able  to  track.  Unfortunately,  a succession 
of  severe  rain  and  wind  storms  had  an  extremely  adverse  effect  on  the 
operation  of  the  COCOA-Cross  telescope  and  not  only  produced  major  inter- 
ruptions in  the  observational  program  and  reduced  its  reliability  during 
synoptic  observation  programs,  but  also  necessitated  the  diversion  of 
financial  and  manpower  resources  away  from  scientific  analysis  toward 
field  maintenance  and  Instrumental  improvements  which  were  made  in  an 
effort  to  Improve  the  reliability  of  the  telescope.  Consequently  the 
TPT  proved  to  be  a valuable  back-up  to  fill  gaps  in  the  COCOA-Cross 
observations. 

Improvements  to  the  COCOA-Cross  have  been  made  both  in  its  field 
hardware  (d.c.  control  wiring,  switches,  general  strengthening)  and  in  its 
pointing  and  data  collection  hardware,  bringing  it  closer  to  our  goal  of 
producing  an  instrument  capable  of  sustained,  uninterrupted  IPS  observa- 
tion (Section  4). 

Synoptic  IPS  observations  were  undertaken  on  both  telescopes 
during  the  summer  of  1976  (17  May-25  August)  and  the  fall  of  1977 
(18  November-13  December).  These  observations,  in  which  measurements 
were  made  of  the  IPS  spectrum,  coupled  with  the  earlier  (1974)  measure- 
ments of  the  IPS  index  (rms  variation  in  intensity),  have  lead  to  the 
conclusion  that,  using  the  techniques  already  developed,  prediction  of 
geomagnetic  activity  is  possible  with  a lead  time  of  about  24  hours 
(Section  2). 

Additional  results  Include  the  association  of  spectral  broad- 
ening with  solar  wind  density  enhancements  near  earth  (Gotvole  et  al.t 
1978 ; and  Appendix  A),  the  detection  of  a solar  wind  velocity  and  tur- 
bulence latitudinal  gradient  out  of  the  ecliptic  plane  (Mitahell,  1978atb ; 
and  Appendix  B),  and  the  extension  of  our  testing  of  the  IPS  prediction 
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technique  to  include  a flare-driven  plasma  disturbance  (from  the  great 
flare  of  November  22,  1977),  and  its  associated  geomagnetic  activity 
(Section  3). 
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Finally  theoretical  considerations  discussed  in  this  report 
| suggest  a new  technique  which  might  extend  the  prediction  lead  time 

significantly,  to  2 or  3 days  (Section  5), 

I 

I 

I 


2. 


SCIENTIFIC  RESULTS 


IC  was  determined,  both  from  attempts  to  Interpret  scintillation 
index  data  in  tersm  of  prediction  (Erakine  at  al.t  1978)  and  from  mathe- 
matical wdellng  (Mitohell  and  Roelof , 1976) , that  in  order  to  explore 
fully  the  feasibility  of  using  IPS  for  prediction  of  geomagnetic  dis- 
turbances, power  spectral  analysis  of  the  intensity  fluctuations  was  highly 
desirable.  The  need  for  spectra  was  based  both  upon  the  dependence 

of  the  frequency  of  the  intensity  fluctuations  on  distance  to  the  scat- 
tering medium,  which  was  essential  information  for  prediction,  and  upon 
the  Increased  ability  to  detect  radio  frequency  Interference  (a  serious 
source  of  contamination  in  the  scintillation  index  data) . The  exploita- 
tion of  this  frequency  dependence  is  discussed  briefly  in  the  Data 
Analysis  Section,  and  more  completely  in  Qotwola  et  al.t  1978. 

Therefore  the  construction  of  a high  speed  data  acquisition 
system  (DAS)  was  completed  at  APL,  resulting  in  a system  which  functioned 
nearly  flawlessly  and  which  included  the  capability  of  simultaneously 
recording  10  sample  per  second  data  both  on  the  COCOA-Cross  and  on  the 
University  of  Maryland  TPT  radio  telescopes.  This  system  is  described 
more  fully  in  Appendix  A.  Also  at  APL  an  extensive  software  development 
effort  was  undertaken  to  reduce  the  raw  data  to  power  spectra  suitable 
for  analysis  (example.  Appendix  A,  Figure  8). 

The  first  data  set  obtained  was  recorded  using  the  APL  DAS  in 
June,  July  and  August  of  1976,  reduced  and  analyzed  at  APL,  and  prediction 
feasibility  possibilities  were  reported  in  the  Interim  report.  Appendix  A. 

A previous  attempt  in  using  IPS  for  prediction  concentrated  on 
the  analysis  of  a 1974  data  set  which  included  only  synoptic  scintillation 
index,  with  no  spectra.  This  was  found  to  be  very  noisy  data,  and  the 
use  of  spectra  in  the  1976  data  was  a marked  improvement  over  the  technique 
of  using  scintillation  index.  In  fact,  an  attempt  was  made  to  use  index 
for  prediction  Independent  of  the  spectral  data  for  the  1976  observations, 
with  noteable  lack  of  success.  In  analyzing  the  1976  spectral  data,  we 
concluded  at  that  time  that  given  continuous  synoptic  observations  on  a 
large  number  of  sources,  particularly  sources  in  the  ecliptic  plane. 


prediction  of  geomagnetic  disturbances  could  be  achieved  with  moderate 
confidence  about  one  day  In  advance.  Since  that  time  we  have  conducted 
another  observation  program  during  October  1977  using  the  only  North- 
South  arm  of  the  COCQA-Cross  on  a reduced  number  20)  radio  sources. 

The  analysis  of  this  data,  which  was  also  recorded  using  the  JHU/APL 
data  acquisition  system,  is  contained  In  this  report.  Section  3,  and 
yields  the  same  conclusion,  however  now  documented  for  a large  solar 
flare  initiated  disturbance  (not  encountered  In  the  1974  or  1976  data 
which  were  dominated  by  corotating  structure).  Again,  a total  fluctua- 
tion power  parameter  similar  to  scintillation  Index  failed  to  produce 
as  positive  a result  as  the  spectral  analysis  -did. 

The  1976  data  set  also  yielded  four  scientific  papers:  Cronyn 
et  al. , 1978;  Gotwola  at  al. , 1978  and  Mitchell,  1978a,b.  The  contents 
of  the  Gotwola  et  al. , 1978  paper  were  similar  to  the  analysis  in 
Appendix  A,  in  which  broadening  of  the  spectra  was  correlated  with  solar 
wind  structure  and  geomagnetic  activity,  while  scintillation  index  was 
found  to  be  so  noisy  as  to  be  of  little  or  no  use  in  prediction. 

To  give  an  idea  of  the  other  scientific  results  which  came  out 
of  the  1976  data  set,  we  describe  the  technique  and  results  from  Mitchell 
1978a  In  Appendix  B.  The  results  from  the  1977  observations  will  now  be 
discussed  in  the  following  section. 
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3.  SCIENTIFIC  ANALYSIS  OF  1977  OBSERVATIONS 

During  the  upgrading  of  the  COOOA-Cross  radio  telescope,  there  was 
a period  between  November  18  (day  322)  and  December  13  (day  347),  1977  over 
which  It  was  possible  to  obtain  nearly  continuous  synoptic  observations  on 
up  to  twenty-two  radio  sources  which  yielded  usable  IPS  spectral  data,  using 
just  the  north-south  arm  of  the  array.  This  configuration  resulted  In  a 
broad  east-west  beam  pattern  with  typical  source  transits  of  ~ 40  minutes. 
Although  this  configuration  necessitated  observations  on  a reduced  number 
of  sources,  the  Increased  time  of  observation  of  a given  source  allowed  us 
to  obtain  spectra  with  excellent  statistical  stability.  We  have  analyzed 
the  data  both  in  terms  of  (a)  the  IPS  power  spectrum  and  (b)  the  total 
integrated  scintillation  power. 

a.  Spectrum  Analysis 

The  spectra  were  parametrized  by  V3,  the  frequency  at  which  the 
power  drops  3-dB  below  the  low  frequency  plateau.  Due  to  Fresnel  diffrac- 
tion effects  and  to  source  angular  diameter  dependence  the  high  frequency 
portion  of  the  spectrum  is  dominated  by  the  near-earth  Interplanetary 
medium  (s  0.  3 AU)  while  the  low  frequency  portion  of  the  spectrum  is  more 
sensitive  to  the  more  distant  medium.  The  source  angular  size, 6,  atten- 
uates intensity  fluctuations  due  to  medium  irregularities  of  angular 
sizes  smaller  than  the  source,  i.e.,  fluctuation  frequencies  above  v “ 

V /TI6L.  The  Fresnel  frequency  v,  " V (CTAL)  the  frequency  below 
1 1 1 

which  the  fluctuation  spectrum  is  flat  for  a power-law  electron  density 
turbulence  spectrum,  also  increases  with  decreasing  distance  along  the 
line  of  sight,  L.  In  these  relations,  \ is  the  radio  frequency  wavelength 
and  V is  the  solar  wind  velocity  projection  perpendicular  to  the  line  of 
sight.  We  attempt  to  exploit  this  distance  dependence  by  using  a measure 
of  the  width  of  the  flat,  low  frequency  portion  of  the  spectrum  (V3)  as 
an  indicator  of  the  distance  between  earth  and  a solar  wind  disturbance 
which  dominates  the  scintillation.  This  technique  is  discussed  in  detail 
in  Gotuola  et  al.t  1978. 
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In  Figure  1 we  compare  the  evolution  of  the  width  of  V3  with  Kp 
index.  The  radio  sources  used  are  listed  from  top  to  bottom  in  order  of 
their  observation  times  from  the  beginning  to  the  end  of  the  UT  day,  and  the 
magnitude  of  V3  is  plotted  logarithmically.  When  two  values  of  V3  are 
plotted  on  a given  day,  the  spectrum  was  what  we  term  "two-component",  l.e., 
it  falls  off  from  the  low  frequency  plateau  approximately  as  a power  law, 
then  flattens  due  to  an  enhancement  at  higher  fluctuation  frequencies  before 
falling  off  again  to  che  white  noise  background  level. 

Ideally  we  would  have  an  ecliptic  point  source  to  the  east  and 
sunward  of  the  earth  with  24  hour/day  monitoring  of  the  scintillations  to 
watch  the  change  in  the  spectrum  caused  by  the  evolution  in  the  density 
fluctuation  profile  associated  with  approaching  solar  wind  disturbances. 
Instead  we  have  finite  diameter  sources  each  of  which  we  can  watch  only 
for  about  40  minutes/day  and  none  of  which  lie  in  the  optimum  position 
(at  this  time  of  year)  for  early  detection  of  a corotating  solar  wind 
disturbance.  In  fact  the  better  scintillators  (3C48,  3C123,  3C144,  3C196, 
3C216)  lie  in  the  anti-sunward  direction  and  to  the  west.  Therefore  in 
exploring  feasibility  of  prediction  we  will  generally  examine  the  width 
of  the  spectra  of  these  western  sources  after  the  onset  of  solar  wind- 
triggered  geomagnetic  activity  to  look  for  the  signature  of  the  solar-wind 
disturbance  as  it  corotates  away  from  the  earth.  Since  at  34  MHz  we  primarily 
respond  to  the  density  increase  preceding  the  rise  in  velocity  in  a solar 
wind  stream-stream  interaction  region  (Erskine  et  al.,  1978),  the  region  is 
not  extended  and  we  expect  similar  signatures  in  the  west  for  a receding 
disturbance  to  that  in  the  east  for  an  approaching  disturbance. 

The  first  event  we  will  examine  is  probably  not  corotating.  A 
flare  on  day  326:09  (November  22)  is  most  likely  responsible  for  the  geo- 
magnetic disturbance  which  began  with  an  SSC  on  day  329.  In  Figure  2, 
we  have  displayed  the  sources  observed  in  ecliptic  coordinates  to  show  the 
detection  of  this  disturbance.  Open  circles  are  narrow  spectra,  closed 
circles  are  broad  spectra,  and  half-filled  circles  are  for  spectra  with 
moderate  high  frequency  enhancement.  A data  gap  for  all  sources  but  3C144 
on  day  329  complicates  the  analysis,  but  we  would  expect  to  see  the  event 
first  in  the  sunward  hemisphere.  Of  the  sources  in  that  direction,  only 


7 


three  (3C273,  3C298  and  3C324)  were  monitored  just  prior  to  the  SSC.  Of 
these,  3C273  does  not  show  a strong  signature,  but  It  does  have  a merglnaliy 
significant  high  frequency  component  on  day  328  (not  recorded  in  Figure  1) , 
consistent  with  the  approach  of  the  disturbance.  3C298  shows  an  Increase  In 
from  day  327  to  328  and  3C324  shows  a relatively  high  on  day  328,  both 
consistent  with  ~ 1 day  prediction  of  the  event.  It  can  be  seen  In  Figure 
2a  that  all  three  of  these  sources  lie  at  elongation  angles  small  enough 
that  strong  scattering  may  have  an  Important  effect  on  the  spectrum.  In 
either  weak  or  strong  scattering  we  would  expect  broadening,  however.  On 
day  329  we  have  data  on  only  one  source,  3C144.  The  observation  took  place 
before  the  SSC,  and  as  3C144  Is  In  the  anti-solar  direction,  the  narrow 
spectrum  observed  is  consistent  with  the  disturbance  not  yet  having  reached 
the  earth.  Figure  2a  shows  the  observations  for  24  hours  before  the  SSC. 

After  the  SSC,  the  Kp  Index,  plotted  In  Figure  1,  rises  and 
remains  elevated  for  half  a day  and  somewhat  enhanced  until  early  on  day 
331  over  which  period  virtually  all  the  sources  observed  display  scintil- 
lation with  broad  spectra  (see  also  Figure  2b) , indicating  that  the  earth 
Is  immersed  In  solar  wind  turbulence.  Although  the  Kp  Index  decays  early 
on  day  331  and  the  Interplanetary  disturbance  has  probably  propagated 
beyond  the  earth,  the  spectra  for  the  anti-solar  sources  3C13,  3C48,  3C68.2, 
3C144,  3C161  and  3C196  remain  broad  on  day  331  (Figure  2c)  as  the  distur- 
bance no  longer  engulfs  the  earth  but  Is  still  near  the  earth  (s  0.3  AU). 

This  argues  for  this  technique's  ability  to  detect  a disturbance  up  to  0.3 
AU  from  earth.  On  day  332  the  spectrum  for  3C48  (a  small  diameter  source 
and  therefore  sensitive  to  turbulence  at  larger  distances  than  most  of  the 
other  sources)  still  exhibits  a high  frequency  enhancement,  but  all  the 
other  sources  return  to  narrow  spectra. 

Figure  2 shows  the  evolution  of  the  response  to  this  disturbance 
and  Indicates  how  we  can  detect  turbulence  in  the  solar  wind  before  it 
arrives  at  earth  (Figure  2a) , when  It  surrounds  the  earth  (Figure  2b) , 
and  after  It  has  gone  past  the  earth  (Figure  2c). 

There  is  a small  enhancement  in  Kp  on  day  334,  which  was  preceded 
by  several  hours  by  a broad  spectnas  on  3C409  looking  to  the  east  (Figure  1). 
Actually,  this  Is  the  second  consecutive  day  that  the  spectrum  for  3C409 
Increased  in  width,  though  It  Is  hard  to  say  whether  all  of  the  broadening 


can  b«  attributed  to  the  disturbance  which  apparently  corotated  past  the 
earth  on  day  334 , since  this  event  Is  closely  followed  by  a such  larger 
enhancement  of  Kp  and  SSC  late  on  day  335  and  continuing  to  the  end  of 
day  336.  Of  the  eastward- looking  sources,  on  which  we  would  expect  to 
see  broadened  spectra  as  a corotating  disturbance  approaches  from  the  east, 
only  3C459  sees  the  disturbance  early  on  day  335,  before  It  reaches  earth. 

We  do  not  have  data  on  the  other  eastern  sources  (4021.53,  3C409,  3C380), 
and  we  cannot  go  more  than  one  day  back  to  look  for  the  signature  of  this 
disturbance  due  to  the  dominance  at  that  time  of  the  day  334  disturbance. 

On  day  336  nearly  all  the  sources  display  broad  spectra,  reflecting 
the  proximity  of  the  turbulence  which  engulfs  the  earth  at  this  time.  The 
next  day  the  spectra  return  to  normal  width  with  the  exception  of  3C216  and 
3C254  which  still  appear  to  be  responding  to  the  solar  wind  disturbance 
which  has  corotated  to  the  west,  again  demonstrating  our  ability  to  detect 
the  disturbance  even  when  it  has  corotated  away  from  the  earth. 

After  this  event  the  data  is  not  very  continuous.  There  is  an 
enhancement  in  Kp  index  on  day  339  which  is  well  associated  with  a broaden- 
ing In  the  spectrum  of  3C298,  and  another  enhancement  In  Kp  on  days  345, 

346,  and  a little  on  347,  which  causes  the  spectra  of  several  sources  to 
broaden.  Of  these  latter  3C216  shows  a broader  spectrum  on  day  347  (after 
the  Kp  Index  has  diminished  somewhat)  than  it  does  on  day  345  at  the  height 
of  the  event,  making  It  a good  candidate  for  prediction  during  the  time  of 
year  It  lies  to  the  east. 

In  an  attempt  to  characterize  the  broadening  of  a spectrum  due  to 
enhancement  at  the  high  frequency  end  in  a more  general  way,  we  devised  In 
Figure  3 a qualitative  display  of  the  relative  change  in  the  width  of  the 
spectrum  for  a given  source  from  one  observation  to  the  next.  In  this 
figure,  we  show  an  Increase  (decrease)  in  spectral  width  from  observation 
1 to  observation  2 by  an  arrow  directed  from  (to)  the  day  of  observation 
1 to  (from)  the  day  of  observation  2.  A double-ended  arrow  means  no  change, 
and  a heavy  arrow  means  a large  change.  Observations  on  a given  day  run 
horizontally;  on  a given  source,  vertically,  and  the  sources  are  listed 
from  left  to  right  in  order  of  their  observation  on  a given  day.  Periods 
of  disturbed  geomagnetic  field  (Kp  > 3)  are  marked  at  the  beginning  by  a 
D,  and  followed  by  a Q (for  quiet,  Kp  £ 3)  at  the  time  when  the  disturbance 
dies  away. 
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This  characterisation  removes  the  confualon  generated  by  one**  end 
two-component  spectra,  and  since  it  is  the  change  in  spectral  width  from 
one  obeorvation  to  the  neat  which  is  recorded  in  thie  approach  to  prediction, 
rather  than  tha  absolute  width  of  the  spectrum,  long  term  trends  such  aa 
elongation  angle  dependence  are  minimised. 

A glance  at  Figure  3 will  tell  which  days  are  of  lntereet.  Any 
day  which  shows  a large  number  of  arrow-heads  pointed  toward  it  would  be 
expected  to  be  an  "all-sky"  day  (a  term  coined  by  Krakina  at  al.t  1978,  to 
describe  a day  when  the  earth  is  immersed  in  turbulence  and  most  sources 
which  can  scintillate  do  scintillate,  in  all  directions).  Thus  on  this 
plot  days  330  and  331,  334,  336,  345  and  347  stand  out  as  "all  sky"  days, 
when  we  would  expect  geomagnetic  activity.  For  prediction,  we  must  examine 
the  figure  more  closely,  being  careful  to  separate  eastern  and  western 
sources  and  looking  for  increasing  width  of  spectra  in  the  east  and 
decreasing  width  in  the  west  for  corotating  disturbances,  or  observing 
separately  the  change  in  width  of  sunward  and  anti-sunward  sources  for 
radially  propagating  disturbances. 

On  day  328  the  sunward  source  3C273  shows  a decrease  in  width 
from  the  day  before,  in  spite  of  the  small  high  frequency  component 
referred  to  previously.  This  decrease  in  width,  which  is  mainly  in  the 
low  frequency  portion  of  the  spectrum,  may  be  due  to  the  presence  of 
more  turbulence  to  the  west  on  day  327,  not  seen  at  earth.  There  was 

1 

only  one  observation  made  on  day  329,  and  that  was  in  the  anti-sunward 
direction  where  3C144  gave  a narrow  spectrum.  A few  hours  later  there 
was  an  SSC,  and  broad  spectra  were  recorded  on  moat  sources  observed 
during  days  330  and  331.  In  addition,  many  of  the  sources  in  the  anti- 
sunward  direction  (3C48,  3C196,  3C216,  3C238)  show  decreasing  spectral 
widths  over  the  next  two  days,  consistent  with  the  increasing  distance 
from  earth  to  the  disturbance. 

There  is  an  all  sky  increase  on  day  334  due  to  the  small  distur- 
bance that  day,  but  any  trend  towards  decreasing  width  of  the  spectra  of 
western  sources  is  obacurred  after  one  day  by  the  arrival  of  a large 
disturbance  causing  an  SSC  late  on  day  335.  This  disturbance  results  in 
an  all  sky  day  on  day  336,  and  several  western  sources  for  irfilch  we  have 
data  (3C48,  3C216,  3C238,  3C254,  3C263. 1)  show  decreasing  spectral  widths 


up  to  day  338.  Of  tha  other  western  sources  for  which  we  do  not  have 
data  continuously  during  this  period,  only  3C123  Is  not  consistent  with 
this  Interpretation. 

Again  early  on  day  345  the  Kp  Index  Is  elevated,  and  the  arrows 

Indicate  an  all-sky  enhancement.  As  we  can  see,  after  the  missing  data  on 

day  346  we  continue  to  show  spectral  broadening  on  day  347.  This  may  be 
detection  of  the  disturbance  which  triggered  the  day  345  event  or  it  may 
be  a separate  disturbance.  The  Kp  Index  remains  somewhat  elevated  through- 
out the  period,  though  it  is  lower  on  day  347  than  on  day  345. 

b.  Scintillation  Power 

Now  we  will  examine  Figure  4,  which  is  a plot  of  where 

Imax  is  the  maximum  I observed  for  a given  source  over  the  course  of  the 

observation  period,  and  I Is  a measure  of  the  normalized  spectral  fluctu- 
ation power: 

J VRMS2-3.15  V23  Hz2 
1 V23  Hz 

2 2 
VRMS  is  the  Integrated  power  over  the  full  spectrum  and  V23  Hz  is  the 

power  Integrated  between  2 Hz  and  the  anti-aliasing  filter  at  3 Hz,  a portion 

of  the  spectrum  containing  virtually  100X  white  noise.  The  3.15  multiplier 

is  an  empirical  constant  chosen  such  that  I ■ 0 for  white  noise.  This 

quantity  has  the  disadvantage  of  being  sensitive  to  the  telescope  response, 

so  that  the  calculated  I's  should  be  normalized  by  the  height  of  the  source 

transit  on  that  day.  This  was  not  done  for  this  data  set  due  to  the  large 

measurement  uncertainty  in  measuring  the  transit  heights  for  many  of  the 

weaker  sources.  We  did  compare  the  daily  response  on  the  strong  source 

3C144  and  normalize  the  3C144  points  to  it.  We  found  a general  decrease 

in  sensitivity  of  the  instrument  over  the  period  by  about  u factor  of  3, 

but  it  was  gradual  and  did  not  seriously  distort  short  term  (3-4  day) 

trends  In  the  data. 
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If  we  look  at  3C196,  3C216,  3C238  and  3C254,  thay  all  allow  a 
decreasing  trend  from  day  330  to  day  332  or  333,  which  My  be  linked  to 
the  scattering  region  moving  farther  from  earth  and  producing  less  total 
fluctuation  in  intensity  due  to  the  Increasing  Influence  of  source  slss 
in  quenching  the  scintillations.  Other  sources  with  I-dependance  consistent 
with  this  trend  are  3C459,  3C48,  3C68.2  and  3C123.  The  exceptions  are 
3C144  and  3C161,  but  a close  examination  of  the  3C144  spectrum  on  day  332, 
the  most  drastic  deviation  from  the  trend,  revealed  a strong  possibility 
of  contamination  by  radio  frequency  Interference.  The  slight  Increase 
between  days  330  and  331  on  3C161  and  3C144  is  probably  due  to  spatial 
structure  of  the  disturbance,  l.e. , the  turbulence  may  have  increased 
between  3C161  and  3C196  on  day  330.  Note  that  unlike  the  spectral  width 
data,  there  is  no  apparent  indication  of  the  approach  of  this  event  in 
the  I-dependence  of  3C298  and  3C324  later  on  day  328.  This  lack  of 
signature  may  be  due  to  strong-scattering,  which  is  a much  more  difficult 
phenomenon  to  recognize  and  Interpret  in  the  total  fluctuation  power  than 
it  is  in  the  spectral  data. 

On  day  334  there  was  an  Increase  in  (normalized)  I on  most  of  the 

sources  which  we  observed  that  day,  consistent  with  the  rise  in  Kp  index  on 

that  day,  and  then  a decrease  on  day  335  leading  into  a large  increase  on 
nearly  all  sources  on  day  336.  The  one  exception  to  this  pattern  for  which 
we  had  data  continuously  from  day  333  to  337  was  3C459.  This  source  dis- 
played increases  on  each  day  from  day  334  to  day  336,  which  we  interpret 
as  first  the  day  334  increase  seen  by  the  other  sources,  but  then  instead 

of  returning  to  a relatively  low  level  as  in  the  case  of  the  western 

sources,  its  fluctuation  power  is  enhanced  further  on  day  335  by  the 
influence  of  the  approaching  corotating  turbulence  which  causes  the  day 
336  enhancement  for  the  other  sources.  Both  the  day  334  and  336  enhance- 
ments are  "all-sky"  days  (as  in  Erakine  et  al.t  1978),  and  the  pattern 
following  the  day  336  enhancement  does  not  show  a mono tonic  decrea8e  as 
clearly  as  the  V3  plot  does. 

The  remainder  of  the  data  is  too  sparse  to  analyze  to  any  degree 
of  precision.  One  event  should  be  mentioned,  though.  Most  of  the  sources 
were  enhanced  on  days  345  and  347  (there  was  no  data  on  day  346),  consistent 
with  the  increase  in  Kp  on  day  345  and  the  continued  Influence  on  the  western 
sources  of  the  associated  solar  wind  disturbance  two  days  later  on  day  347. 
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Wc  conclude  that  In  looking  at  the  data  set  In  terns  of  both 
spectral  shape  (v^)  and  scintillation  power  (I)  as  nsasured  by  the  OOCQA- 
Cross  in  1977,  Interplanetary  scintillations  at  34  MHs  were  useful  for 
prediction  of  geomagnetic  activity  at  beat  about  one  day  in  advance, 
provided  the  radio  source  locations  are  optimum  for  prediction,  l.e., 
sunward  and  to  the  east,  but  not  in  strong  scattering.  It  maybe  that 
instrumental  and  analytical  techniques  not  employed  here,  such  as  that 
discussed  in  Section  5,  will  improve  the  lead  time  of  the  prediction 
using  34  Wi  radio  scintillations. 


ft  » 
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a.  Field  Hardware 

For  orientation  purposes,  simplified  schematics  of  the  north  (or 
south)  and  east  (or  west)  telescope  arms  are  shown  in  Figures  5 and  6.  It 
should  be  understood  that  the  COCOA-Cross  was  originally  constructed  with 
a grant  of  $30,000  in  1972  from  NASA,  and  numerous  economies  in  materials 
and  construction  techniques  had  to  be  employed  which  resulted  in  poor 
reliability  and  vulnerability  to  long  term  deterioration.  For  example, 
all  the  d.c.  lines  which  controlled  the  telescope  phasing  diode  switches 
were  grouped  in  surplus  multiconductor  cables  which  lay  directly  on  the 
ground.  After  four  years  of  exposure  to  strong  sunlight,  high  temperatures, 
rain  and  submersion  in  alkali  mud,  the  cable  insulation  was  cracked.  Con- 
sequently extensive,  dif ficult-to- locate  shorts  developed  after  each  rain. 
The  diode  switches,  consisting  of  a diode,  3 capacitors  and  2 resistors, 
were  constructed  on  small  acrylic  plastic  blocks  and  dipped  in  epoxy. 

Again,  after  four  years  of  severe  environmental  exposure  the  epoxy  coating 
had  been  largely  sand-blasted  away,  causing  the  diodes  to  crack  and  com- 
ponent connections  to  corrode  and  fall. 

Thus,  on  the  bacis  of  the  condition  of  the  d.c.  control  lines 
ana  switches,  and  the  over-riding  requirements  for  a high  reliability 
Instrument  to  take  24  hour/day,  365  day/year  IPS  observations,  we  decided 
to  install  new  d.c.  control  lines  and  design  and  install  new  phasing  control 
diode  switches.  The  new  d.c.  control  lines  are  open  wire  lines,  entirely 
elevated  above  the  ground  (see  Figure  7).  These  lines  are  not  only  above 
any  reasonably  predicted  water  level  but  also,  because  the  lines  are  open, 
shorts  and  breaks  can  be  readily  located  and  fixed  without  having  to  cut 
into  heavy  jacketed  multiconductor  cable  (as  was  necessary  for  the  old  d.c. 
control  lines).  Even  if  the  lines  should  be  submersed  in  an  extraordinary 
flood,  they  can  dry  out  quickly,  unlike  the  multi-conductor  cables  which 
tended  to  trap  water  inside  the  outer  jacket.  A total  of  approximately  30 
miles  of  new,  elevated,  open  wire  d.c.  control  line  was  installed  by  August 
1977. 


All  of  the  diode  phasing  switches  were  also  replaced  by  a total 
of  approximately  1600  opaque,  fully  encapsulated  switches  (Figure  8). 


i m 


In  addition  to  the  switches  and  control  lines,  we  also  Installed 
preamplifiers  at  the  output  of  each  bank  of  8 antenna  elements  to  improve 
the  overall  signal-to-noise  ratio  of  the  telescope  (Figure  9).  Measurements 
taken  on  just  the  north-south  area  in  November-December  1977  showed  that  the 
signal-to-noise  ratio  in  the  IPS  power  spectrum  was  3 db  higher  than  for  the 
entire  telescope  before  the  new  bank  and  buffer  amplifiers  were  added.  The 
signal-to-noise  ratio  for  the  entire  telescope  is  6 db  higher  than  for  the 
north-south  area  alone  so  we  anticipate  that  the  total  improvement  in  signal- 
to-noise  ratio  will  be  about  9 db.  A total  of  16  preamplifiers  were  installed 
in  the  north-south  arm,  and  32  in  the  east-west  arm.  Circuitry  to  enable  on- 
off  control  of  each  individual  preamplifier  was  also  installed. 

b.  Control  Hardware 

In  the  original  d.c.  control  line  system  each  type  of  switch  (for 
example,  all  north-south  arm  line/fan  position  2 switches)  was  connected  to 
its  own  indivudual  d.c.  control  line  which  was  either  turned  on  to  forward  bias 
the  switches,  or  turned  off  which,  for  the  silicon  diode  switches,  was  effect- 
ively back-biasing  the  switch.  That  is,  the  control  voltage  was  unipolar. 

However,  to  reduce  the  number  of  new  d.c.  control  lines  by  a 
factor  of  2 and  thereby  reduce  installation  and  maintenance  time  and  expense, 
we  decided  to  employ  a bipolar  switch  control  voltage  on  each  line  for  which 
one  polarity  would  serve  to  forward-bias  one  member  (i.e.  N-S  line/fan  2)  of 
a pair  of  switches  and  back-bias  the  complementary  member  (i.e.  N-S  line/fan 
1);  reversing  the  polarity  reversed  the  switch  pair  combination  (i.e.  to  N-S 
line/fan  2 off,  N-S  line/fan  1 on).  The  old  beam-steering  control  system  was 
designed  to  provide  only  unipolar  switching  voltages.  In  an  effort  to  provide 
a simple  unipolar-to-bipolar  switching  voltage  conversion  we  attempted  to  use 
a simple  reed-relay  interface.  The  solenoids  were  controlled  by  the  original 
unipolar  switching  voltage  while  the  armature  switched  between  + and  - 48 
volts.  Unfortunately  the  reed  relays  had  a poor  reliability  record  (due  in 
part  to  the  highly  capacitive  d.c.  control  line  main  feed  cable),  and  very 
occasional  short  circuits  (i.e.  once  during  the  night  in  a wind  storm)  would 
blow  the  protective  fuses,  rendering  the  beam  steering  system  useless.  Although 
the  observations  of  November-December  1977  were  taken  using  the  reed  relay  inter- 
face it  necessitated  an  inordinate  amount  of  attention  and  d.c.  control  line 
current  monitoring  to  insure  proper  operation. 
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We  therefore  felt  that  it  was  essential  to  undertake  the  design 
of  new  electronic  solid-state  switches  which  would: 

1.  Provide  protection  against  external  shorts  via  current 
limiting,  not  fusing,  so  that  occasional  shorts  would  have  no  permanent 
(i.e.,  fuse-blowing)  effect  on  the  beam  steering  system; 

2.  Provide  immediate  visual  and  electrical  flags  to  indicate 
deviations  in  total  switch  current  draw  (both  excesses,  or  shorts,  and 
deficiencies,  or  breaks); 

3.  Be  easily  inter faceable  to  the  Poly  8080-based  processor. 

The  new  switch  drivers  not  only  have  features  (1.)  - (3.)  but 
also  have  a controllable  disable  to  turn  off  1 or  more  selected  d.c. 
control  lines  for  test  purposes,  and  a built-in  dual  beam  capability  to 
provide  for  essentially  Instantaneous  switching  between  2 different  beams 
without  having  to  transfer  complete  beam  pointing  information  each  time. 
One  of  the  switch  driver  cards  is  shown  in  Figure  10. 

The  switch  driver  cards  are  controlled  by  the  Poly  8080-based 
microprocessor  system  through  a single  master  interface  card  (Figure  11) 
and  a switch  driver  control  card  (1  for  each  6 switch  drivers;  see 
Figure  12). 


- 16  - 


[ 

[ 

[ 

[ 

[ 

[ 


c.  Poly  8080-Based  Microprocessor 

Repeated  mechanical  failures  of  the  card  reader  and  paper  tape 
punch  and  their  associated  continuing  maintenance  costs  and  the  major 
interruptions  to  the  observing  program  caused  by  d.c.  control  line  faults 
and  the  necessity  of  carrying  out  long,  tedious  diagnostic  procedures  on 
the  control  lines  compelled  us  to  consider  the  use  of  a completely  non- 
mechanical  telescope  control  and  data  acquisition  system  based  upon  the 
new  inexpensive  microprocessor  chips.  The  system  requirements  were  four- 
fold: 

(1)  Replace  the  card  reader  and  directly  control  the  upgraded 
main  beam  switching  network; 

(2)  Sense  malfunctions  in  the  telescope  itself; 

(3)  Replace  the  punch  unit  and  record  the  data  on  a more 
reliable  audio  cassette  recorder; 

(4)  Upline,  or  transfer,  to  the  NOVA  computer  system  the  data 
stored  on  the  audio  cassette. 

The  final  hardware  configuration  of  the  new  system  consists  of  a 
Poly  8080-based  micro  computer  with  16  K bytes  of  RAM,  a cassette  modem, 
a video  monitor,  a cassette  recorder,  a keyboard,  and  a parallel-serial 
I/O  handler.  The  micro  system  is  connected  via  a 4800  band  serial  line 
to  the  NOVA  computer  system.  An  8-level  vectored  priority  interrupt  and 
a real-time  clock  are  part  of  the  hardware  architecture. 

An  overall  block  diagram  of  the  beam-pointing  and  data  acquisition 
system  is  shown  in  Figure  13. 

The  software  system  was  developed  in  assembly  language  to  optimize 
memory  use  and  to  take  full  advantage  of  the  interrupt  hardware.  The  custom- 
ized real-time  operating  system  used  a software  sidereal  clock  to  control  the 
various  functions  of  the  telescope,  to  gather  data  and  record  it,  and  to 
respond  to  external  commands  to  perform  other  tasks  - such  as  receive  a new 
observing  program,  replay  the  data  accumulated  on  the  audio  tape,  reset  the 
clocks,  modify  the  beam  positioning,  page  through  the  observing  list,  etc. 

The  microcomputer  is  not  a slave  to  the  NOVA;  rather,  the  two  systems  com- 
municate as  peers.  This  allows  the  telescope  to  function  independently  of 
the  data  processing  and  scientific  calculations  being  carried  out  off-line 
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by  Che  NOVA.  The  micro  is  a system  dedicated  to  telescope  monitoring  and 
control,  and  data  acquisition,  while  the  NOVA  serves  as  a general  purpose 
computational  facility. 

The  observing  procedure  is  as  follows:  Sources  are  selected 
from  the  Cambridge  IPS  catalog  or  other  lists  of  sources  which  are  known  to 
be  of  small  angular  size  and  display  IPS  at  34.3  MHz.  The  coordinates  are 
entered  into  a NOVA  program  called  ALLIN1  which  produces  processed  direct 
coordinates,  elliptic  and  galactic  coordinates,  and  the  beam  steering 
switch  settings.  These  data  are  recorded  on  a magnetic  tape  for  transfer 
at  a convenient  time  to  the  memory  of  the  microcomputer  control  system. 

The  telescope  observes  the  sources  on  the  observing  schedule,  which  may 
include  up  to  200  sources.  The  radio  intensity  and  phase  information 
channels  are  digitized  and  the  data  temporarily  stored  in  the  computer 
memory.  At  the  present  data  rate  of  one  channel/5  sec.,  the  memory 
buffer  is  damped  onto  audio  cassette  tapes  approximately  once  an  hour. 

A standard  30  minute  cassette  tape  provides  a data  storage  capacity  of 
at  least  5 full  days.  At  a convenient  time  this  audio  cassette  data 
tape  is  rewound  and  read  back  into  the  microcomputer  which  uplines  the 
data  via  the  serial  connector  to  the  NOVA  computer.  This  process  does 
not  interrupt  the  normal  observing  of  the  telescope.  The  data  is  rec- 
orded in  NOVA  compatible  format  for  off-line  analysis.  During  the  ob- 
serving procedure  at  regular  Intervals  the  different  control  lines  are 
checked  for  undercurrent,  caused  by  defective  array  switches  or  broken 
lines,  and  overcurrent,  often  caused  by  short  circuits.  This  error 
information  is  presented  on  a CRT  monitor  to  Inform  the  observer  of 
possible  maintenance  problems. 

The  on-site  data  reduction  procedures  available  to  the  ob- 
server have  been  improved  significantly.  Not  only  has  the  inconvenience 
of  paper  tape  handling  been  eliminated,  but  a video  terminal  Installed. 

Data  can  be  quickly  examined,  manipulated,  and  the  final  results  displayed 
on  a Versatec  printer/plotter . New  programs  have  been  prepared  to  aid  in 
the  Fourier  analyis  of  the  data  and  in  theoretical  calculations  pertinent 
to  the  solar  wind  and  IPS.  This  capability  formerly  was  not  present  at 
the  site  itself. 
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AND 

We  emphasize  that  we  have  by  no  means  exploited  the  full  perfor- 
mance and  operational  capability  of  the  telescope,  and  that  with  a longer 
data  base  of  higher  signal- to-noise  radio  observations  on  more  sources, 
which  the  telescope  Is  now  capable  of,  the  lead  time  may  be  extended. 

That  is,  we  feel  that  the  feasibility  of  making  predictions  based  on  IPS 
activity  has  not  been  fully  explored  during  this  contract  period,  but 
the  physical  reliability  of  the  telescope  has  been  drastically  improved, 
array  malfunction  monitoring  instrumentation  has  been  installed,  and  data 
acquisition  and  processing  is  sufficiently  highly  automated  that  the  small 
staff  of  scientific  personnel  stationed  at  Clark  Lake  on  the  COCOA-Cross 
project  can  realistically  expect  to  keep  up  with  the  large  number  of  daily 
observations  and  parameters  (4  data  parameters/sources ; ~ 125  sources/day 
going  to  ~ 200/ day). 

All  of  the  new  d.c.  control  lines,  amplifiers,  diode  phasing 
switches  and  switch  drivers,  the  Poly  8080-based  microprocessor  system 
and  the  control,  data  acquisition  and  data  processing  software  were 
developed  and  installed  only  because  they  were  absolutely  essential  to 
the  project  mission  of  studying  the  feasibility  of  using  IPS  to  forecast 
geomagnetic  disturbances.  The  size  of  the  telescope,  harshness  of 
environment,  volume  of  observational  data  and  lack  of  skilled  manpower 
made  necessary  the  hardware  and  software  additions  and  improvements. 

It  was  unfortunate  that  the  initiation  of  the  feasibility  study  coin- 
cided with  the  wettest  and  most  destructive  weather  of  the  past  16  years 
(see  Appendix  C) . 

We  have  previously  (Roelof  et  al .,  1977-AFGL  Fiiul  Rc['ort) 
discussed  the  fact  that  for  a source  of  angular  size  8,  at  a given 
frequency  v in  the  power  spectrum,  contributions  to  fluctuation  power 
will  be  strongly  attenuated  at  distances  > L: 

L - V /H0v  - 0.18  AU(V  / A00  km  s_1)/(P"v) 

X X 

so  that  for  typical  V / A00  1,  6"  *>  1 and  v ■ 0.5  Hz,  L^x  0.35  Al'. 

The  measurements  made  to  date,  i.e.,  scintillation  index  for  the  chart 
record  data  and  frequency  scale  for  the  power  spectra,  are  strongly 
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influenced  by  near-by  turbulent  structure.  However,  If  absolute  power 

2 

spectral  density  normalised  by  (source  response)  is  measured  at  low 
frequencies,  l.e.,  ~ 0.15-0.2  Hz,  then  the  maximum  distance  of  response 
Is  extended  to  ~ 0.7  AU  for  6 ~ 1.0"  which  means  that  the  lead  time  for 
corotating  density  enhancement  warnings  is  extended  to  ~ 3 days.  As  we 
resume  operation  with  the  up-graded  COOOA-Cross  facility  we  shall  be 
monitoring  total  fluctuation  power  at  ~ 0.2  Hz  to  exploit  this  possi- 
bility of  extending  lead  time. 
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6. 


SUMiARY 


The  program  undertaken  to  study  the  feasibility  of  prediction 
using  IPS  was  conducted  as  a collaboration  between  the  University  of  Iowa 
group  at  the  COCQA-Cross  Radio  Telescope  installation  at  Clark  Lake  Radio 
Observatory  which  provided  manpower  for  up-grading  the  COOOA-Cross  array 
and  the  use  of  the  instrument  observations,  and  The  Johns  Hopkins  Univer- 
sity/Applied Physics  Laboratory  group.  The  latter  assumed  responsibility 
for  data  collection  and  analysis,  including  mathematical  modelling,  con- 
struction of  a data  acquisition  system,  observers  on  site,  data  reduction 
and  analysis,  and  publication  of  scientific  results  both  individually  and 
Jointly  with  members  of  the  Iowa  group. 

Damage  to  the  COCOA  Cross  as  a result  of  abnormal  rainfall, 
and  the  continued  vulnerability  of  the  telescope  to  rain  damage  accom- 
panied by  extended  periods  of  down  time,  made  it  necessary  to  replace 
the  original  diode  phasing  switches  and  multiconductor  d.c.  control 
lines  with  opaque,  fully  encapsulated  switches  and  elevated  open  wire, 
d.c.  control  lines.  These  changes,  combined  with  the  limitations  and 
failures  of  the  card-reader  controlled  main  beam  steering  system  and 
the  punched  paper  tape  data  acquisition  system,  compelled  us  to  develop 
and  produce  new,  bipolar,  solid-state  switch  drivers  controlled  by  a 
microprocessor  which  also  serves  to  acquire  the  low  rate  data  previously 
taken  on  the  punched  tape  system. 

Bank  and  buffer  amplifiers  which  have  significantly  improved 
the  signal-to-noise  ratio  of  the  telescope  have  been  installed. 

The  improvements  in  field  hardware  and  control,  monitoring 
and  data  acquisition  hardware  and  software  now  make  it  possible  to 
undertake  a far  more  intensive  prediction  feasibility  study  than  was 
possible  at  the  beginning  of  the  contract  period.  The  telescope  is 
now  also  far  less  vulnerable  to  rain  and  wind  damage,  and  far  easier 
to  maintain. 

As  a result  of  the  damage  to  the  telescope  caused  by  severe 
weather  conditions,  the  extensive  data  set  anticipated  at  the  inception 
of  this  program  was  never  realized.  Two  separate  sets  of  observations 
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ware  obtained,  in  the  summer  of  1976  and  the  fall  of  1977.  Aa  continuous 
synoptic  observation  Is  essential  to  the  application  of  113  to  prediction, 
the  usefullness  of  the  observations  obtained  was  further  limited  by  data 
gaps  caused  by  telescope  malfunctions  and  radio  frequency  Interference 
(both  of  which  should  be  greatly  reduced  after  the  recent  improvements). 

The  feasibility  of  using  IPS  for  prediction  is  therefore  based  upon 
rather  undefinitive  data,  and  it  is  our  recommendation  that  further 
observations  be  carried  out  to  confirm  the  results  of  this  study,  namely 
that  prediction  using  this  technique  Is  feasible  under  favorable  conditions 
with  lead  times  up  to  and  sometimes  exceeding  24  hours.  During  the  second 
observation  period  we  extended  the  application  of  the  technique  beyond 
purely  co-rotating  solar  wind  disturbances  to  include  flare  associated 
disturbances  which  resulted  in  geomagnetic  activity. 

The  use  of  the  spectra  of  the  Intensity  fluctuations  has  been 
fundamental  to  the  understanding  of  this  technique.  Furthermore,  appli- 
cation of  the  spectral  broadening  technique,  developed  in  the  course  of 
this  study,  has  led  to  significant  scientific  results  including:  the 
association  of  the  major  component  of  the  IPS  index  response  and  spectral 
broadening  at  34  MHz  with  density  enhancements  in  the  solar  wind;  and  the 
detection  of  a positive  heliolatitude  gradient  in  solar  wind  velocity 
associated  with  a possible  negative  gradient  in  solar  wind  turbulence. 

Though  it  is  desireable  to  continue  the  acquisition  and  reduc- 
tion of  some  data  for  spectral  analysis,  the  work  (both  analytical  and 
theoretical)  conducted  during  this  study  has  provided  us  with  the  under- 
standing and  insight  into  the  frequency  dependence  of  the  intensity 
fluctuations  necessary  to  make  informed  decisions  on  measuring  more 
easily  obtained  parameters  which  will  be  useful  for  prediction.  There- 
fore on  the  basis  of  observations  to  date,  and  supporting  theoretical 
analvsis,  we  have  decided  to  record  not  only  the  total  source  response, 
fluctuation  power  from  0.1  to  1.5  Hz  and  first  spectral  moment  (which 
like  1b  a measure  of  spectral  width) , but  also  aource-response- 
normalized  fluctuaticu  power  at  relatively  low  fluctuation  frequencies, 
l.e.,  ~ 0.1-0. 3 Hz.  We  anticipate  that  such  measurements  will  yield 
predictive  measurements  with  IPS  activity  leading  geomagnetic  activity 
by  upwards  of  two  days. 
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FIGURE  CAPTIONS 

FIGURE  1 Plot  of  spectrum  3 dB  frequency  on  a logarithmic  scale  versus 
day  of  the  year  for  twenty  radio  sources.  Kp  Index  Is  plotted 
at  the  bottom  for  the  same  period,  versus  day  of  the  month. 
Sudden  storm  commencements  (SSC's)  are  Indicated  by  a triangle 
both  on  the  Kp  plot  and  on  the  log  (V3)  plot  between  the 
sources  whose  observations  bracket  the  SSC.  Dashed  lines 
connect  points  separated  by  more  than  one  day.  Points  in 
parentheses  are  the  second  component  of  possible  two-component 
spectra. 


FIGURE  2 Ecliptic  longitude  plots  for  the  day  prior  to  and  two  days 

after  an  SSC  on  day  329:1200.  Filled  circles  indicate  large 
increases  in  V3  or  broad  spectra  (for  that  source),  half-filled 
indicate  moderate  increases  In  v-j  or  intermediate  width  spectra, 
and  open  circles  indicate  decreases  in  or  narrow  spectra. 

On  the  day  prior  to  the  SSC  (a)  only  sunward  sources  show 
broadening  in  response  to  the  disturbance  inside  1 AU. 

Following  a data  gap,  day  330  and  early  day  331  (b)  shows  an 
all-sky  day,  with  broad  spectra  in  all  directions  as  the  earth 
is  immersed  in  turbulence,  and  the  following  day  (c)  the  broad 
spectra  are  almost  entirely  early  in  the  day  and  in  the  anti- 
sunward  direction,  responding  to  the  disturbance  as  it  propa- 
gates beyond  the  earth  into  deep  space. 

FIGURE  3 Plot  of  the  (qualitative)  change  in  width  of  the  spectrum  from 
one  observation  to  the  next  for  19  radio  sources  over  the 
November-December  1977  observation  period.  Light  arrows 
indicate  a moderate  change  in  spectral  width,  pointing  toward 
the  day  with  the  broader  spectrum.  Heavy  arrows  Indicate  a 
large  change  in  spectral  width,  and  double-ended  arrows  indicate 
little  or  no  change.  D and  Q designate  the  beginning  of  periods 
of  disturbed  and  quiet  geomagnetic  field,  respectively.  Sudden 
storm  commencements  are  indicated  by  a filled  triangle  placed  on 
the  day  of  occurrence,  between  the  source  observations  which  fall 
just  prior  to  and  just  after  the  SSC. 
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FIGURE  4 Linear  plot  of  normalized  fluctuation  power  (1/1,^^)  versus 

time,  with  Kp  index  plotted  below.  Dashed  lines  connect  points 
separated  by  more  than  one  day.  Sudden  storm  commencements  are 
indicated  by  filled  triangles  placed  on  the  day  of  occurrence 

between  the  source  observations  Immediately  preceeding  and 

following  (in  time)  the  SSC.  This  is  a measure  Blmllar  to 
scintillation  index,  but  is  not  normalized  by  the  source 
intensity. 

FIGURE  5 Simplified  schematic  of  East  (or  West)  arm  branch  feed  system. 

All  transmission  lines  are  balanced  open  wire  line.  Lines  are 
Identified  by  number:  within  the  column  north-south  lines  1-3; 
between  columns,  east-west  lines  1-5.  Antenna  elements  (V)  are 
approximately  33  m long  and  oriented  east-west.  North-south 
separation,  6.55  m.  Phase  switching  system  same  as  in  north- 
south  arm.  Two  control  lines  are  required  for  each  of  the  3 
north-south  lines,  3 are  currently  required  to  accomplish  east- 
west  pointing,  and  2 (not  shown)  are  required  to  accomplish 
north-south/east-west  arm  phasing  for  a total  of  11  control 
lines.  The  original  unipolar  switches  were  replaced  by  998 
fully  encapsulated  switches  controlled  by  bipolar  switching 
voltages.  Bank  and  buffer  amplifiers  were  also  installed.  The 
switch  and  amplifier  tasks  were  made  possible  by  AFOSR  support. 

FIGURE  6 Simplified  schematic  of  North  (or  South)  arm  branch  feed  system. 

All  transmission  lines  are  balanced  open  wire  line.  Lines  are 
identified  by  number  (1-8)  on  left.  Antenna  elements  (V)  are 
approximately  63  m long  and  oriented  east-west.  North-south 
element  separation  is  6.55  m.  Appropriate  phase  for  reversing 
switch  (0°  or  180°)  is  determined  by  polarity  of  line  reversing 
switch  control  voltage,  Vc,  and  similarly  for  the  phase  of  the 
fan  switch  (0°  or  90°).  Each  type  of  switch  (for  example  all 
32  north  and  south  arm  line  2,  fan  switch  2 diode  switches)  is 
controlled  by  its  own  individual  d.c.  control  line.  Since  fan 
and  reversing  switch  control  voltages  are  required  for  each  of 
lines  1-6,  3 real-time  delay  tap  voltages  are  required  on  line 
7 and  5 on  line  8,  and  3 control  lines  are  required  to  set  the 
phasing  switches  (not  shown)  between  the  north  and  south  arms, 
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FIGURE  7 


a total  of  23  d.c.  control  lines  plus  a common  return  line  are 
required  for  north-south  arm  beam  pointing.  The  original  uni- 
polar switches  were  replaced  by  608  fully  encapsulated  switches 
controlled  by  bipolar  switching  voltages.  Bank  and  buffer 
amplifiers  were  also  installed.  Both  switch  and  amplifier 
tasks  were  made  possible  by  AFOSR  support. 

New  d.c.  control  lines  reflecting  sun  light.  The  lines  are 
elevated  at  least  16  inches  off  the  ground. 


FIGURE  8 New,  fully  encapsulated  phasing  switches.  Switch  cable  is 
connected  to  new  d.c.  control  line. 

FIGURE  9 Bank  amplifier  at  output  of  East  arm  bank. 

FIGURE  10  Switch  driver  board,  4 drivers/board. 

FIGURE  11  Master  interface  board  to  connect  Poly  88  microprocessor  to 
control  boards. 

FIGURE  12  Control  board  for  switch  drivers.  Each  board  contains  6 

switches  (1^  switch  drive  boards).  It  is  connected  through 
a bidirectional  data  line  to  a master  interface.  Each  of 
the  13  control  boards  has  a unique  address  set  on  the  binary 
slide  switches. 

FIGURE  13  Block  diagram  of  beam-pointing  and  data  acquisition  system.  Source 
ID  and  coordinates  are  entered  into  the  NOVA  which  then  completes 
current  celestial,  ecliptic  and  galactic  coordinates;  phase  up  and 
transit  times;  switch  settings;  and  miscellaneous  control  settings 
are  computed  by  a NOVA  program  called  ALLIN1.  A printer  listing 
and  cassette  tape  containing  data  downlined  to  the  Poly  88  are 
generated  by  the  NOVA.  Switch  settings  are  distributed  by  the 
Poly  88  through  an  interface  board  to  control  boards  and  the 
switch  driver  boards.  If  phasing  switch  current  draw  is  either 
lighter  or  heavier  than  present  thresholds  a flag  is  set  which 
is  passed  to  the  Poly  88.  The  flags,  time  and  4 channels  of 
data  are  recorded  on  a digital  cassette  tape  for  convenient 
upllnlng  to  the  NOVA  for  final  processing  and  plotting. 
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This  interim  scientific  report  covers  the  first  year  (1  April  1976  - 
31  March  1977)  of  the  research  conducted  under  AFOSR  contract  NPP-75-157 
"Prediction  of  Solar  Particle  Events  and  Geomagnetic  Activity  Using  Inter- 
planetary Scintillation  Observations  from  the  Iowa  COCOA-Cross  Radio  Telescope". 

In  September  1975  a two  year  program  of  observations  of  interplanetary 
radio  wave  scintillations  aimed  at  the  prediction  of  solar  particle  events  and 
geomagnetic  activity  was  proposed  to  APOSR.  This  proposal  called  for  the  instal- 
lation at  the  Clark  Lake  Radio  Observatory  of  a digital  data  acquisition  system 
which  would  be  designed  and  built  at  APL.  Further  it  was  proposed  to  evaluate 
the  feasibility  of  the  use  of  the  observations  for  prediction  of  solar-terres- 
trial disturbances. 

Observations  were  taken  May-August  1976  at  CLRO  and  reduced  and  ana- 
lyzed at  the  Applied  Physics  Laboratory.  Interplanetary  scintillation  power 
spectra  were  examined,  and  a significant  new  prediction  signature  for  approach- 
ing solar  wind  turbulence  was  identified.  This  signature,  which  characterizes 
the  shape  of  the  spectrum,  contains  more  information  than  the  commonly-used 
scintillation  index,  which  contains  the  band-passed  total  power  of  the  spectrum. 
We  Intend  to  pursue  the  study  of  spectral  signatures  as  a diagnostic  of  plasma 
turbulence  approaching  the  earth.  During  the  summer  of  1976,  the  turbulence 
detected  by  spectral  analysis  was  well-associated  with  geomagnetic  disturbances 
caused  by  turbulent  solar  wind  streams  (measured  by  spacecraft)  impinging  on 
the  magnetosphere.  The  turbulence  detection  preceded  the  geomagnetic  distur- 
bance by  at  least  one  day. 

2.  SUMMARY  OF  OPERATIONS:  APRIL  1976  - MAY  1977 

Our  proposal  in  September  1975  called  for  the  Installation  at  the 
Clark  Lake  Radio  Observatory  of  a digital  data  acquisition  system  within  6 
months  after  the  commencement  of  funding.  In  fact  the  data  acquisition  was 
installed  in  only  1.5  months.  This  compression  in  schedule  was  possible  because 
of  the  provision  by  APL  of  independent  research  and  development  funds  which 
allowed  construction  of  the  data  acquisition  system  to  begin  before  funds  were 
received  from  AFOSR. 
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Observations  were  conducted  continuously  on  the  COCOA-Cross  Radio 
Telescope  (34.3  MHz)  from  mid  Kay  1976  until  late  August  1976.  Of  the  ~ 150 
radio  sources  observed  each  day,  a subset  of  41  sources  proved  favorable 
enough  in  location,  flux  and  angular  diameter  to  produce  measureable  scintil- 
lation power  (see  Table  1) . The  sources  are  shown  in  ecliptic  coordinates  in 
Figure  1.  Of  these,  15  consistently  rendered  spectra  with  sufficient  signal- 
to-noise  ratios  to  be  used  in  the  analysis  of  the  scattering  in  the  inter- 
planetary medium.  Observations  on  the  COCOA-Cross  were  supplemented  with 
simultaneous  observations  at  38  MHz  on  the  University  of  Maryland  TPT  array, 
also  at  CLRO,  using  the  same  data  acquisition  system.  At  the  conclusion  of 
the  observation  period  it  was  clear  that  the  data  set  could  be  substantially 
improved  by  upgrading  the  array,  in  order  to  improve  the  signal  to  noise  ratio, 
reduce  suseptlbility  to  radio  frequency  interference,  and  increase  the  reli- 
ability of  the  system  during  adverse  weather  conditions.  This  upgrading  is 
discussed  more  fully  in  the  next  section. 

Interpretation  of  the  data  began  in  the  fall  of  1976  and  is  contin- 
uing at  present.  A brief  discussion  of  our  preliminary  results  is  presented 
in  Section  5 and  a list  of  publications  in  Section  7. 

3.  INSTRUMENTATION 

3.1  Digital  Data  Acquisition  System 

A block  diagram  of  the  data  acquisition  system  constructed  at  APL 
for  use  in  this  study  is  shown  in  Figure  2.  This  system  allows  up  to  4 analog 
signals  to  be  simultaneously  digitized  and  recorded  on  magnetic  tape.  A number 
of  optional  sample  rates  are  available  up  to  the  maximum  rate  of  10  samples/sec. 
In  addition,  each  analog  channel  has  an  event  channel  associated  with  it  which 
allows  a single  bit  of  information  to  be  recorded,  such  as  marking  the  instant 
when  the  antenna  was  phased  up  for  a new  source.  All  timing  information  is 
derived  from  the  highly  accurate  (crystal  controlled)  serial  time  code  which 
is  supplied  from  the  University  of  Maryland's  time  code  generator. 

Since  power  spectra  are  being  calculated  from  the  10  SPS  data  collected 
by  this  system,  it  is  particularly  important  to  minimize  the  amount  of  aliasing 
that  occurs  above  the  5 Hz  Nyqulst  frequency.  This  is  accomplished  by  placing  a 
low  pass  filter  in  front  of  the  analog-to-dlgltal  converter.  This  filter  is  a 
4 pole  Chebyshev  design  with  a cutoff  frequency  of  2.9  Hz,  and  an  attenuation  at 
the  Nyquist  frequency  of  20  dB.  This  perhaps  overly  conservative  design  virtu- 
ally guarantees  freedom  from  spectral  ambiguities  due  to  aliasing. 
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TABLE  1 — Sources 

for  which  scintillation 

index  was  routinely  i 

I 

3C13 

3C147 

3C286 

3C33 

3C154 

3C287 

I 

* 3C48 

3C161 

3C293 

3C55 

* 3C186 

3C295 

r 

* 3C68.2 

* 3C196 

* 3C298 

i * 

3C79 

* 3C216 

3C310 

r 

3C82 

3C225 

* 3C380 

I 

1. 

3C89 

3C252 

* 4C+21.53 

I* 

3C103 

* 3C254 

3C409 

L 

3C111 

* 3C263.1 

3C436 

1 1' 

3C123 

3C268.4 

3C438 

3C125 

* 3C270.1 

3C456 

3C134 

* 3C273 

* 3C459 

I r 

* 3C144 

* 3C280 

* Source  for  which  the  power  spectrum  was  frequently  calculated. 
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3C161 


Fig.  1 Subset  of  41  sources  observed  by  the  COCOA-Cross  radio  telescope. 

Sources  plotted  with  an  open  circle  consistently  yielded  useful  spectra, 


liagram  of  the  data  acquisition  system  used  to  obtain  data 
le  COCOA-Cross  and  TPT  arrays. 


Despite  the  almost  dally  occurrence  of  blowing  sand,  the  COCOA- 
Cross  equipment  van  is  not  yet  equipped  with  any  particular  environmental 
controls.  Since  this  would  lead  to  a greatly  shortened  life  for  the  digi- 
tal tape  recorder,  permission  was  obtained  from  the  University  of  Maryland 
to  place  the  bulk  of  the  data  acquisition  system  In  the  more  controlled 
environment  of  the  Maryland  observatory.  The  "front  end"  (filters  and  A/D 
converters)  was  then  Installed  In  the  COCOA-Cross  van  and  a digital  link 
established  between  the  two  parts  of  the  system.  Although  not  Indicated 
on  the  block  diagram.  It  is  possible  to  run  either  all  channels  remote 
(COCOA-Cross) , all  channels  local  (University  of  Maryland) , or  two  channels 
remote  and  two  channels  local.  This  flexibility  allows  for  the  possibility 
of  a coordinated  observing  program  between  the  COCOA-Cross  and  the  University 
of  Maryland  TPT  variable  frequency  array  when  the  TPT  is  available  for  IPS 
studies. 


3.2  COCOA-Cross  Array  Upgrading 

There  has  been  a continuing  problem  with  long  down  times  of  the 
COCOA-Cross  following  the  occasional  heavy  rain  storms  which  have  occured  at 
Clark  Lake.  The  problems  have  been  confined  to  the  DC  control  line  cables 
used  for  controlling  the  pointing  of  the  telescope,  and  the  radio  frequency 
phasing  control  diode  switches.  The  roughly  100  conductor-miles  of  multi- 
conductor cable  were  Installed  at  ground  level  and  were  therefore  totally 
submerged  under  highly  corrosive  and  conductive  alkali  water  and  mud  when 
the  lake  was  flooded.  The  diode  switches  (512  in  the  N-S  arm  and  896  in  the 
E-W  arm)  were  constructed  on  plexiglass  plates  and  given  a thin  coating  of 
epoxy  which  has  been  severely  sand-blasted,  resulting  in  exposure  and  damage 
to  the  components  of  some  switches,  and  the  generation  of  sporadic  noise, 
particularly  after  a rain.  The  old  switches  and  control  lines  are  shown  in 
Figures  3 and  4 and  the  Improved  switches  and  lines  in  Figures  5 and  6. 

As  discussed  elsewhere  in  this  report  a productive  synoptic  IPS 
spectral  measurement  program  was  carried  out  from  May  to  August  1976.  Obser- 
vations had  been  suspended  up  until  the  summer  because  of  damage  from  heavy 
rains  in  March.  However,  in  early  September  the  program  had  to  be  suspended 
again  after  the  "dry  lake"  was  flooded  by  heavy  rains  which  accompanied  Hurri- 
cane Kathleen.  In  early  October,  before  the  lake  hail  dried  out  enough  to 


View  of  E arm  column  looking 
North.  Note  DC  control  line 
cables  on  ground  at  base  of 
pole. 


Figure  4 Excellent  picture  of  a pair 
of  old-8tyle  switches.  Note 
hew  visible  components  are 
(very  thin  coating  of  epoxy) 
and  exposed  individual  switch 
wires . 


resume  observations,  the  Lake  was  flooded  for  a third  time  in  1976  by  heavy 
rains.  The  resulting  damage  to  the  DC  control  lines  was  so  heavy  and  exten- 
sive that  we  decided  to  concentrate  all  on-site  COCOA-Cross  staff  operations 
on  replacement  of  the  control  lines  and  switches. 

All  the  multiconductor  DC  control  lines  have  now  been  replaced  by  a 
network  of  copperweld  open-wire  lines  which  are  at  least  12  inches  above  the 
ground,  and  all  of  the  1408  diode  switches  have  been  replaced  by  switches  which 
are  totally  encapsulated  in  epoxy  (see  Figures  5 and  6).  Final  checkout  of  the 
reconfigured  array  was  delayed  somewhat  due  to  damage  incurred  in  early  March 
1977  as  the  result  of  a very  severe  wind  storm  with  gusts  reaching  125  mph.  At 
this  writing  this  damage  has  been  corrected  and  a full  checkout  of  the  new  con- 
trol lines  and  switches  is  underway.  In  addition,  the  physical  construction  of 
the  array  has  been  strengthened  in  many  areas  so  as  to  significantly  improve  its 
resistance  to  wind  damage. 

It  is  expected  that  synoptic  IPS  observations  will  begin  in  June  and 
will  no  longer  be  interrupted  by  extensive  down  times  following  heavy  rains 
and  very  high  winds. 

4.  DATA  REDUCTION 

All  raw  data  are  transcribed  from  the  7 track  tape  written  at  Clark 
Lake  to  9 track  tape  for  study  at  APL.  The  main  beam  response  and  scintillation 
power  in  the  0.1  to  1.5  Hz  band  are  then  digitally  calculated  for  each  observa- 
tion and  printed  on  a Versatek  printer-plotter  (Figure  7).  The  scintillation 
index  is  then  calculated  as  described  in  Figure  7.  These  computer  generated 
plots  have  the  great  advantage  over  our  former  (prior  to  May  1976)  strip-chart 
type  of  presentation  in  that  the  computer  can  first  scan  the  entire  observation 
and  select  the  best  one  of  the  available  standard  dynamic  ranges  for  displaying 
each  channel.  In  this  way,  the  scaling  errors  of  our  old  scheme,  which  occurred 
when  the  width  of  the  track  became  comparable  to  the  deflection,  have  been  essen- 
tially eliminated. 

Examination  of  the  time  history  plots  mentioned  above  allows  rapid  iden- 
tification of  sources  which  are  likely  to  yield  useful  power  spectra.  The  power 
spectrum  analysis  consists  of  first  taking  out  long-term  trends  by  passing  the 
data  through  a digital  high  pass  filter  (4  pole,  matched  Z Butterworth  design) 


Figure  5 New  switch.  Components  cere  totally  encapsulated,  in  opaque  epoxy. 
Switch  leads  are  protected  inside  cable  jaaket. 


Figure  6 A new  switch  assembly  (description  as  per  Figure  S ).  DC  control 
lines  are  elevated  far  above  any  possible  water  level  and  are 
open  so  that  there  are  no  cable  joints  which  require  sealing. 
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Fig.  7 Scintillation  power  squared  (SP2)  and  main  beam  (MB)  response  for 
three  successive  transits  of  3C48  on  1976,  day  191.  Unbarred 
variables  are  2.5  second  averages  which  are  plotted  chiefly  to  aid  in 
the  identification  of  interference.  Barred  variables  are  50  second 
triangle  weighted  averages  which  are  scaled  by  hand  (in  minor  chart 
divisions)  and  then  combined  with  the  scale  factor  (KM)  just  to  the 
left  of  each  SP2  curve  to  yield  the  scintillation  index.  For  example, 
for  the  leftmost  set  of  curves  the  scintillation  index  is:  m = 1.789  * 
(v/45.0)/27.0  = 0.44.  The  Q values  shown  are  quality  factors  which 
are  helpful  when  summarizing  a large  number  of  observations. 
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with  a cutoff  frequency  of  0.02  Hz.  These  hlghpass  data  are  then  broken  into 
256  point  segments  (25.6  sec)  and  a fast  Fourier  transform  is  performed  on  each 
segment.  The  average  magnitude  squared  of  each  Fourier  component  is  then  cal- 
culated and  corresponding  components  from  each  256  point  transform  are  added 
together  to  give  the  desired  spectral  estimates.  On  those  sources  which  are 
observed  more  than  once  per  day  all  of  the  spectral  estimates  are  further  com- 
bined into  one  composite  spectral  estimate  with  Improved  stability  (Figure  8). 

In  order  to  further  increase  the  statistical  stability  of  the  spectral  estimates, 
the  spectra  are  smoothed  with  a Hanning  spectral  window.  This  technique  is 
equivalent  to  a combination  of  Bartlett  and  Hanning  smoothing.  For  example, 
for  the  spectrum  shown  in  Figure  8,  which  represents  9 minutes  of  duta,  there 
are  112  degrees  of  freedom. 

Two  curves  are  plotted  in  Figure  8.  The  upper  curve  is  the  raw 
uncorrected  spectrum  which  consists  of  a scintillation  component  plus  a white 
noise  component  due  to  system  noise.  The  sharp  cut-off  above  3 Hz  is  due  to 
an  anti-aliasing  filter.  Since  in  our  observations  the  scintillation  component 
of  the  spectrum  has  always  been  found  to  be  negligible  above  2 Hz,  we  estimate 
the  white  noise  amplitude  by  averaging  across  the  2 to  3 Hz  range  and  subtract 
that  estimate  from  each  of  the  spectral  points.  This  yields  the  lower  curve 
which  for  convenience  has  been  renormalized.  The  low  frequency  droop  below 
about  0.04  Hz  is  due  to  the  digital  high-pass  filter. 

Three  empirical  parameters  are  tabulated  from  each  of  the  many  spectra 
examined.  These  are:  and  v^,  the  frequency  at  which  the  corrected  spectrum 

falls  3 and  6 dB  respectively  below  the  low  frequency  plateau,  and  a,  the  high 
frequency  slope  of  the  power  law.  The  frequencies  and  are  not  taken 
directly  from  the  spectral  points  but  instead  are  read  from  an  eye  estimate 
of  the  best  fitting  curve  to  the  data  points. 


5.  DATA  ANALYSIS 

The  reduced  data  were  analyzed  using  comparisons  with  models  of  inter- 
planetary scintillation  and  upplied  to  a study  of  the  feasibility  of  using  them 
for  the  prediction  of  possible  geomagnetic  activity  resulting  from  the  inter- 
action of  the  earth's  magnetosphere  with  solar  wind  disturbances.  Displaying 
the  spectra  on  a log-log  plot  (example.  Figure  8),  they  typically  displayed  a 
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Scintillation  frequency  (Hz) 


Fig.  8 Power  spectrum  of  the  scintillating  radio  source  3C48  on  1976,  day 
191. 
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flat  plateau  from  very  low  frequencies  to  a turn-over  into  a nominal  power- 
law  slope,  the  turn-over  occurring  at  anywhere  from  less  than  0.1  Hz  up  to 
0.5  Hz.  Both  the  frequency  and  the  rate  at  which  the  spectrum  turns  over 
into  the  power  law  provide  information  about  the  velocity  of  the  scattering 
medium  and  its  density  turbulence  profile.  The  slope  of  the  power-law  also 
provides  information  about  the  turbulence  profile.  Therefore,  the  spectra  are 
parametrized  by  recording  the  power-law  slope  ( a ) , the  frequency  at  which  it 
drops  to  3 dB  below  the  plateau  (v^) , and  the  frequency  at  which  it  drops  to 
6 dB  below  the  plateau  (v^) . These  three  parameters  then  serve  as  useful 
characterizations  of  the  spectra  for  the  analysis  of  large  amounts  of  data. 

A plot  of  versus  time  (Figure  9)  for  the  15  sources  on  which  we 
had  reasonable  continuity  demonstrates  the  utility  of  this  sort  of  analysis 
for  the  prediction  of  the  arrival  at  earth  of  solar  wind  disturbances.  The 
days  on  which  solar  wind  stream-stream  interaction  regions  arrived  at  the 
earth,  as  determined  from  IMP-7  and  8 spacecraft  data.  Figure  10,  are  indi- 
cated in  Figure  9,  and  shaded  on  the  figure  about  those  sources  which  were 
located  favorably  for  prediction.  The  most  useful  source  during  this  obser- 
vation period  for  the  purpose  of  predicting  the  arrival  of  disturbances  was 
3C273.  It  lies  to  the  east  of  the  sun,  from  which  direction  corotating  dis- 
turbances come, and  it  is  at  a large  enough  elongation  angle  so  that  it  is  not 
in  strong  scattering.  It  also  lies  virtually  in  the  ecliptic  plane  so  that 
disturbances  confined  to  low  solar  latitudes  lie  in  the  line  of  sight,  unlike 
the  case  with  higher  latitude  sources  such  as  3C263.1  and  3C270.1. 

The  two  events  for  which  we  have  the  best  coverage  are  day  182  and 
day  197.  In  each  of  these  cases,  the  spectra  on  3C273  displayed  increases 
in  (broadening  in  the  spectra)  at  least  1 day  before  the  arrival  of  the 
disturbance  at  earth.  This  increase  in  is  what  one  would  expect  on  the 
basis  of  modelling  of  the  scattering  in  the  medium,  as  a turbulent  region 
moves  closer  to  the  earth.  A close-up  study  of  the  geometry  and  the  power 
spectrum  for  3C273  leading  up  to  the  geomagnetic  activity  on  day  197  is  shown 
in  Figure  11.  In  addition,  the  day  197  event  appears  to  have  been  seen  by 
other  eastern  sources  (3C270.1,  3C298) , Indicating  that  the  turbulent  region 
extended  above  the  ecliptic  plane  in  the  north.  Unfortunately,  3C216  and 
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Fig.  9 p3  versus  time  for  15  sources  - known  solar  wind  streams,  SC's  and 
abrupt  increases  in  Kp  index  are  indicated.  The  symbols  are  repeated 
(shaded)  to  point  out  the  areas  in  the  data  consistent  with  prediction 
of  events.  On  those  spectra  displaying  two  plateaus  both  values  of  P3 
are  plotted. 
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Fig.  10  Solar  wind  probe  data  ( Lazarus , Solar-Geophysical  Data,  1977),  days 
1 53-244,  showing  recurring  corotating  solar  wind  streams.  Shaded 
areas  correspond  to  known  streams  with  sufficient  IPS  data  coverage 
to  be  used  in  the  prediction  study.  Triangles  are  geomagnetic 
sudden  commencements;  arrows  are  abrupt  increases  in  Kp  index. 
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Fig.  1 1 (a)  Smoothed  scintillation  spectra  for  3C273  on  days  195, 196  and 
197  illustrating  progressive  broadening  as  a corotating  region 
approaches  the  earth. 

(b)  Line  of  sight  to  3C273  projected  onto  the  ecliptic  plane  on  1976, 
day  197.  Also  shown  is  the  approximate  shape  of  the  center  of  a 
corotating  density  enhencement  for  days  195, 196  and  197  at  the 
time  when  3C273  was  observed. 
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3C196  are  both  at  elongations  small  enough  to  be  in  strong  scattering,  for 
which  the  present  prediction  analysis  is  not  well  suited. 

Another  aspect  of  the  data  analysis  is  shown  in  Figure  12.  These 
are  scatter  plots  of  versus  (measured  on  the  same  spectrum),  which  give 
an  Indication  of  the  rate  at  which  the  spectrum  turns  over  into  a power-law. 

The  degree  to  which  the  points  follow  a straight  line  indicates  the  constancy 
of  the  ratio  v^/v^,  and  the  plots  show  that  to  a good  approximation  the  spectra 
represent  universal  curves  whose  locus  is  shifted  back  and  forth  in  frequency 
by  conditions  in  the  scattering  medium. 

Figure  13  is  a display  of  scintillation  index  (m)  versus  time  for 
several  eastern  sources,  and  there  is  no  apparent  trend  in  the  data  (including 
that  for  3C273)  which  could  be  exploited  for  the  purpose  of  prediction.  It  is 
possible  that  this  data  might  hold  more  value  for  prediction  in  future  obser- 
vations with  the  improvements  currently  being  made  on  the  instrument. 

6.  PRELIMINARY  EVALUATION  AND  FUTURE  PLANS  * 

We  have  demonstrated  that  spectral  analysis  of  interplanetary  scin- 
tillations enhances  their  utility  as  a probe  of  the  interplanetary  medium  in 
predicting  the  arrival  of  solar  wind  disturbances  at  earth,  provided  that  many 
sources  are  monitored  to  the  east  of  the  sun.  High-speed  data  collection  and 
reduction  to  spectra  is  therefore  a highly  desirable  step  in  the  process  of 
prediction.  Our  greatest  obstacle  in  the  evaluation  of  this  spectral  technique 
has  been  data  gaps  and  uncertainties  in  the  data  due  to  variable  array  response. 
With  the  improvements  and  upgrading  the  array  is  currently  undergoing  we  anti- 
cipate a more  definitive  evaluation  ensuing  from  another  period  of  data  collec- 
tion. Another  way  in  which  we  hope  to  improve  the  statistics  and  reliability 
of  the  data  is  by  exploiting  more  fully  the  3-beam  east-west  steering  capabil- 
ities of  the  array,  in  order  to  spend  more  time  per  day  on  a given  radio  source. 
This  has  already  been  shown  to  significantly  improve  the  statistics  on  a given 
day's  observation  of  a source  during  the  summer  of  1976  observing  period.  It 
not  only  increases  the  total  amount  of  data  going  into  a spectrum,  but  decreases 
the  probability  of  missing  a source  entirely  due  to  RFI. 
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Fig.  12  Scatter  plots  of  p3  versus  i>g  for  3C48(a)  and  3C273(b).  The  high 
correlation  between  thses  two  parameters  implies  a consistently 
similar  shape  for  the  spectra,  i.e.,  a universal  curve  which  shifts  on 
the  frequency  axis  in  response  to  solar  wind  velocity  and  turbulence 
profile.  The  closer  grouping  of  the  3C273  points  may  be  due  in 
part  to  large  source  angular  diameter  effects.  The  points  in 
parenthesis  for  3C273  are  taken  from  2 component  spectra,  which 
result  from  two  widely  separated  scattering  regions  in  the  inter- 
planetary medium. 
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Fig.  13  Scintillation  index  versus  time  for  8 eastern  sources.  The  shaded 
regions  indicate  the  density  enhancements  proceeding  several  high 
speed  solar  wind  streams.  Abrupt  increases  in  Kp  and  geomagnetic 
sudden  commencements  are  also  indicated. 
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Based  on  the  results  we've  obtained  to  date,  we  have  decided  to 
begin  synoptic  observations  again  in  June  (1977)  with  the  upgraded  array. 
Included  in  the  upgrading  (besides  that  discussed  in  Section  3)  will  be 
high-pass  filtering  of  the  RF  signal  to  reduce  RFI  suseptibility , and  a better 
monitoring  system  to  identify  RFI.  We  intend  to  use  the  3-beam  per  source 
feature,  particularly  in  the  east.  The  data  will  be  reduced  and  analyzed  at 
APL,  and  parallel  theoretical  work  will  be  continued  to  enhance  our  understan- 
ding in  the  observations. 
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APPENDIX  B 

ADDITIONAL  SCIENTIFIC  RESULTS  FROM  1976  OBSERVATIONS 
(Condensed  from  Mitchell,  1978a, b) 

In  an  earlier  paper  (Mitchell  and  Roelof , 1976)  we  presented 
some  calculations  of  Intensity  spectra,  based  upon  simple  assumptions. 
The  electron  density  fluctuation  spectrum  was  assumed  to  have  the  form 


F(k)  - FQ(ko2  + k2)'q/2exp(-k 2/kt2) 


It  was  assumed  that  the  outer  scale  < 2TT / ko ) was  very  much  larger  than 
the  first  Fresnel  zone  radius  L,  where  X is  observing  wavelength 
and  L the  distance  from  the  observer  to  the  scattering  medium),  and 
that  the  inner  scale  ( 2Tl/ k^)  was  near  zero.  Thus  the  inner  and  outer 
scales  have  negligible  effects  on  the  intensity  scintillation  spectrum, 
except  for  the  dependence  of  F on  k for  q > 3,  or  on  k for  q < 3. 

° _4° 

It  was  further  assumed  that  F « r , where  r is  heliocentric  radius. 

0 2 - -2 

This  is  reasonably  close,  but  not  identical  to  assuming  < 6N  >2  « r 
for  an  extended  medium. 

In  order  to  model  the  scattering  from  an  extended  region,  we 
first  divide  the  medium  into  twenty  slabs  along  the  line  of  sight,  each 
of  which  contributes  roughly  equally  to  the  scintillation  index  of  a 
point  source.  Thin-screen  spectra  are  calculated  for  the  middle  of  each 
slab,  and  then  summed  discretely  through  the  medium  to  produce  the 
extended  medium  spectrum.  Besides  the  entire  medium  case  (twenty  slabs), 
other  cases  discussed  in  this  report  include  models  where  the  scatter- 
ing occurs  only  in  the  first  seven  slabs  closest  to  the  earth,  the  first 
thirteen  slabs,  or  a middle  range  (seven  to  thirteen) . These  approxi- 
mations are  more  realistic  than  a single  thin  screen  model  in  allowing 
some  distribution  to  the  turbulence. 

Figure  B1  addresses  the  additional  effect  of  source  size  on  the 
distance  dependence  of  the  scintillation  response.  The  figure  illustrates 
the  relative  contribution  (dm/dL)  to  the  scintillation  index  (m)  as  a 
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function  of  distance  (L)  along  the  line  of  sight  normalized  to  its  maximum 
value,  for  different  elongations  and  source  sizes.  The  locations  of  slabs 
7,  13,  19  and  20  are  indicated.  It  can  be  seen  that  for  large  angular 
diameters,  the  distant  medium  contributes  little  or  nothing  to  the  scin- 
tillation, thus  limiting  the  distance  to  which  the  scintillation  is 
sensitive  as  a probe  of  the  medium.  In  modeling  the  scattering  using 
less  than  the  entire  medium,  the  contribution  to  the  index  will  go  as 
these  curves  do  within  the  slabs  used  and  will  equal  zero  outside  of 
these  slabs.  This  approximation  is  quite  reasonable  when  one  considers 
that  the  density  (and  presumably  the  turbulence)  enhancement  associated 
with  stream-stream  interaction  regions  is  routinely  only  an  order  of 
magnitude  above  quiescent  levels,  and  the  scattering  power  goes  as 
< 6N2  >. 

There  are  several  parameters  which  characterize  the  medium  and 
the  radio  source.  The  solar  wind  velocity  (V)  is  assumed  constant  and 
radial.  The  electron  density  fluctuations  are  assumed  to  have  a power- 
law  spatial  spectrum.  See  Mitchell  and  Roelof , 1976  for  a complete 
discussion  of  the  computations  based  on  this  spectrum  and  of  this 
technique.  The  medium  Irregularities  are  assumed  isotropic,  and  the 
source  is  modelled  as  having  an  isotropic  Gaussian  brightness  distribu- 
tion with  an  e-folding  diameter  6. 

We  further  ignore  receiver  bandwidth  effects  and  limit  the 
treatment  to  an  observing  wavelength  of  8.75  meters.  This  corresponds 
to  the  observing  frequency  (34.3  MHz,  500  KHz  bandwidth)  of  the  Univer- 
sity of  Iowa  COCOA-Cross  radio  telescope. 

The  most  easily  obtained  spectral  parameters  are  the  apparent 
high  frequency  slope  (a)  and  the  turnover  frequency  (v^)  defined  as  the 
frequency  at  which  the  spectrum  drops  to  3 dfi  below  the  Fresnel  plateau. 
The  Fresnel  plateau  is  the  flat  portion  of  the  spectrum  at  low  frequencies 
and  is  defined  operationally  by  drawing  a straight  line  through  the  low 
frequency  points  between  ~ 0.08  Hz  and  the  point  at  which  a smooth  curve 
drawn  through  the  spectrum  falls  off  mono toni cal ly.  A more  complete 
description  of  the  experimental  techniques  employed  are  available  in 
Gotwols  et  al.,  1978. 
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As  th«  errors  inherent  in  a are  greater  than  those  in  and  the 
need  for  accuracy  for  physical  Interpretation  greater  in  a than  in 
(a  has  a smaller  relative  dynamic  range  than  v^)>  we  will  be  primarily 
concerned  with  the  direct  comparison  of  with  the  model  and  will  use 
a only  as  a rough  indicator  of  consistency  between  the  model  and  the  data. 

For  small  sources  (0  £ 2")  whose  spectrum  is  dominated  by  Fresnel 
turnover,  is  by  the  nature  of  its  definition  biased  toward  the  more 
distant  medium.  This  is  because  any  contribution  from  the  distant  medium 
Just  adds  power  at  low  fluctuation  frequencies. 

However,  for  larger  source  sizes  (at  3")  the  scattering  from  the 
most  distant  medium  is  severely  attenuated,  and  so  the  scattering  region  is 
in  effect  more  localized.  In  this  case  becomes  a better  indicator  of 
the  distance  to  the  scattering  medium.  In  the  limit  of  very  large  sources 
(a  5"),  the  dominant  scattering  region  is  near  earth  (£  0.3  AU) . In  this 
case  the  cut-off  of  the  spectrum  is  source  size  dominated. 

We  display  ranges  of  from  theoretical  spectra  at  a range  of 
elongations  (60°  & c £ 120°)  in  Figure  B2,  with  source  size  as  a parameter. 
It  should  be  emphasized  that  for  each  realization  of  the  medium  employed 
here  (with  the  exception  of  the  entire  medium  case),  a change  in  elongation 
angle  of  twenty  degrees  also  represents  a sizeable  change  in  the  shape  of 
the  density  profile,  for  a fixed  choice  of  slab  numbers.  For  example, 
slabs  1-7  are  of  different  thicknesses  and  distances  at  different  elonga- 
tions. In  these  Idealized  cases,  the  variation  in  elongation  can  be  looked 
upon  as  a change  in  the  medium  profile.  In  fact,  for  these  large  elonga- 
tion angles,  the  daily  variability  In  the  turbulence  profile  dominates 
real  elongation  trends  in  the  profile.  Thus  in  comparing  an  observation 
at  a given  elongation  with  a calculated  for  the  proper  source  size, 
one  would  not  necessarily  expect  to  find  a fit  to  the  theoretical  at 
the  same  elongation,  unless  one  closely  approximated  the  actual  turbulence 
and  velocity  profiles.  However,  even  with  the  limited  resolution  of  the 
ranges  of  calculated  in  Figure  B2,  one  can  roughly  infer  source  size 
and  the  configuration  of  the  medium  by  comparing  a sufficient  sampling  of 
experimental  data  with  these  ranges.  Note  that  for  each  realization  of 
the  medium  there  is  a decreasing  progression  of  with  increasing  source 
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six*.  Thu*  two  sources  of  different  sixes  observed  et  the  sane  time  et  the 
same  elongation  would  be  seperated  by  this  dependence.  By  comparing  model 
points  against  synoptic  experimental  data,  one  might  then  be  able  to  deter- 
mine a source's  site  and  the  configuration  of  the  scattering  medium.  One 
obstacle  to  this  approach  is  the  need  to  know  both  the  velocity  of  the 
medium,  and  its  turbulence  profile  along  the  lln*-of-alght.  However,  a 
rough  knowledge  of  source  slxe  should  be  obtainable  even  without  velocity 
Information,  given  enough  synoptic  data,  since  the  deta  should  cluster 
about  mean  turbulence  and  velocity  profiles.  Also,  day-to-day  changes  in 
the  location  of  a traveling  turbulent  structure  should  be  readily  observable 
by  monitoring  the  change  in  width  of  the  spectrum  on  the  same  source  each 
day,  as  the  source  size  is  then  a constant  and  the  velocity  of  the  structure 
should  also  be  nearly  a constant. 

Discussion  - Experiment 

We  present  in  Figure  B3  the  slopes  (or)  and  turnover  frequencies 
(v3)  for  data  taken  during  the  months  of  June  and  July  1976  on  the  sources 
3C48  and  3C273.  The  confidence  in  the  estimate  of  a is  generally  lower 
than  that  for  due  to  the  relatively  lower  confidence  level  of  the  high- 

frequency  points  as  they  approach  the  noise  background  level;  therefore 
the  following  analysis  relies  more  heavily  on  fitting  than  on  fitting  o. 
The  data  in  Figure  B3  were  taken  from  spectra  obtained  at  34.3  MHz  with  the 
University  of  Iowa  COCOA-Cross  radio  telescope.  Additional  observations 
supplementing  the  COCOA-Cross  data  were  taken  at  38  MHz  on  the  University 
of  Maryland  TPT  radio  telescope,  also  at  Clark  Lake  Radio  Observatory. 

It  is  immediately  apparent  that  the  data  from  the  two  sources  lie  in 
distinctly  different  areas  of  the  (a,  v^)  plane.  As  the  scattering  regions 
in  the  Interplanetary  medium  are  physically  separated  for  the  two  sources, 
and  the  sources  have  different  angular  diameters,  this  separation  is  not 
unexpected.  The  3C48  points  ere  sufficiently  high  in  that  regardless 
of  source  size,  we  must  assume  (a)  that  the  medium  is  almost  always 
described  as  having  a relative  excess  of  turbulence  near  the  earth  (con- 
centrated ~ 0.1  AU  in  distance)  or  (b)  that  the  assumed  velocity  of  400  km/ 
sec  is  systematically  exceeded  In  the  medium.  This  is  not  true  of  the  3C273 
points,  which  Implies  a difference  in  the  character  of  the  medium  between 
the  lines  of  sight  to  the  two  sources.  In  Figure  B4,  a subset  of  the  data 
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from  Figure  B3  for  which  we  obtained  corresponding  IPS  velocity  data  from 
UCSD  (John  Harmon , Private  Carrrunication)  is  replotted,  with  the  turnover 
frequency  corrected  to  a standard  of  V ■ 400  km/sec  for  all  observa- 

tions, e.g.,  v3(V)  - v3(400  km/ a)  (V/400).  Included  in  Figure  B4  are 
those  ranges  from  Figure  B2  which  best  fit  each  of  the  sources'  points. 

Note  that  while  the  7-13  slab  ranges  fit  in  v^,  the  or's  calculated  are 
much  higher  than  those  measured  for  angular  diameters  (6)  2 1.0".  It 
is  clear  that  to  a large  degree  the  scatter  of  the  data  from  the  two 
sources  in  in  Figure  B3  (replotted  on  Figure  B4  for  comparison  in 
this  format)  is  attributable  to  the  differences  in  the  velocities  of 
the  scattering  medium  along  the  lines  of  sight  to  the  two  sources. 

For  the  period  during  which  the  observations  overlap  (day  182- 
day  204)  the  average  velocity  for  3C48  is  nearly  200  km/sec  higher  than 
that  for  3C273.  The  corotation  delay  between  the  two  sources'  scattering 
regions  is  nominally  about  3 days  during  this  23  day  overlap  period 
(since  their  ecliptic  longitude  difference  is  ~ 45°)  and  the  velocities 
for  each  source  do  not  fluctuate  by  more  than  ~ 100  km/sec  (i.e.,  no  one 
day  weights  the  average  significantly,  nor  does  the  velocity  for  3C273 
ever  exceed  that  for  3C48)  even  though  there  were  well-developed  corota- 
ting streams  observed  at  this  time. 

If  one  calculates  the  latitudes  of  the  closest  point  of  approach 
to  the  sun  of  the  lines  of  sight  to  the  two  sources,  the  average  gradient 
(assuming  radial  flow)  over  this  23  day  period  is  > 20  km  s xdeg  . 

There  remains  in  Figure  B4  a clear  separation  between  the  locus 
of  the  two  sources'  v^'s.  In  general,  the  3C273  points  lie  at  lower  cut- 
off frequencies  than  the  3C48  points.  Part  of  this  separation  is  due  to 
a difference  in  the  angular  extent  of  the  two  sources.  However,  the  points 
for  3C48  lie  at  higher  turnover  frequencies  than  can  be  explained  by  the 
model  on  the  basis  of  a spherically  symmetric  interplanetary  medium. 
Normally  the  distant  medium  would  add  power  to  the  spectrum  preferentially 
at  low  frequencies,  which  would  result  in  low  (~  0.2  Hz)  turnover  frequen- 
cies. The  higher  frequency  deviation  of  the  3C48  points  may  be  due  to  two 
latitude  effects.  First,  as  the  velocity  data  indicates,  there  is  a strong 
latitudinal  gradient  in  the  solar  wind  velocity.  As  discussed  earlier, 
the  frequencies  at  which  turbulence  at  a given  distance  along  the  line  of 
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sight  contributes  to  the  observed  spectrum  ere  proportional  to  the  bulk 
velocity  of  the  turbulent  region.  If  this  velocity  gradient  oontinuee  to 
higher  latitudes  (>  10s),  it  could  shift  the  fluctuation  power  contributed 
by  the  more  distant  medium  (L  > 1 AU)  to  higher  frequencies;  our  computa- 
tional model  assumed  a constant  radial  velocity.  Second,  the  medium  at 
ecliptic  latitudes  > 10-15°  may  be  less  turbulent  and/or  have  lower  den- 
sities than  the  medium  in  the  ecliptic.  Low  turbulence  levels  would 
produce  less  scattering,  and  if  the  higher  latitudes  have  lower  turbu- 
lence levels,  then  the  turbulence  level  along  the  line  of  sight  to  a non- 
ecliptic source  would  decrease  with  distance  as  the  latitude  of  the  point 
along  the  line  of  sight  increased.  Because  the  distant  medium  contributes 
proportionately  more  power  at  low  frequencies  than  the  near-earth  medium, 
the  low  turbulence  levels  in  the  distant  medium  would  relatively  reduce 
the  low  frequency  power  in  the  observed  spectrum  and  would  increase. 

This  would  explain  the  apparent  absence  of  power  In  the  3C48  spectrum  at 
low  frequencies. 

There  is  another  way  of  seeing  that  both  of  these  effects  contri- 
bute to  the  deviation  of  the  3C48  results  from  the  spherically  symmetric 
model.  If  we  compare  the  distribution  of  the  3C48  points  with  the  calcu- 
lated v^'s  in  Figure  B2,  we  can  see  that  they  come  closest  to  fitting  the 
points  for  which  the  medium  was  modelled  as  extending  only  to  the  13th  of 
the  twenty  screens  of  equal  scattering  power.  They  are  consistent  with  a 
source  size  of  about  ~ 0. 5 to  3 arc  seconds  for  q ■ 2.8.  3C273,  on  the 
other  hand,  is  consistent  with  a source  size  of  — 5 arc  seconds,  and  a 
spherically  synmetrlc  medium.  As  its  ecliptic  latitude  is  4°,  its  scat- 
tering region  is  essentially  in  the  the  sun's  equatorial  plane  and  the 
distribution  of  turbulence  (and  scattering)  along  the  line  of  sight  would 
be  relatively  insensitive  to  latitude  dependence  effects.  Most  of  the 
3C273  points  lie  at  low  turnover  frequencies;  however,  they  fall  over  a 
broad  range  of  slopes  (1.3  a a a 3. 5) . The  implication  here  is  that  the 
thickness  of  the  scattering  medium  is  varying,  the  higher  slopes  resulting 
from  days  on  which  the  scattering  was  dominated  by  a relatively  localized 
disturbance.  Some  of  the  range  in  a is  also  due  to  uncertainties  in  the 
estimate  of  the  slope. 
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Modelling  Che  electron  density  spectrum  es  a power-law  with 
slope  q,  it  was  found  that  the  best  fits  to  the  data  (taken  at  34.3  J®z) 
were  obtained  assuming  q ■ 2.8  (rather  than  q ■ 3.6;  q ■ 3.2  may  also 
produce  a good  fit),  with  source  sizes  of  ~ 0.5  to  3.0  arc-seconds  for 
3C48  and  ~ 5.0  arc-seconds  for  3C273.  It  was  also  found  that  the  spectral 
turnovers  for  3C48  lie  at  higher  3 dB  frequencies  than  can  be  explained 
with  a spherically  symmetric  medium.  We  interpret  this  to  be  the  result 
of  a latitudinal  gradient  In  velocity,  and  when  normalized  by  IPS  velocity 
data,  a continuation  of  that  gradient  In  either  turbulence,  (with  lower 
turbulence  at  higher  latitudes)  or  in  velocity  (with  higher  velocities  at 
higher  latitudes)  or  both. 
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FIGURE  CAPTIONS  FOR  APPENDIX  B 


FIGURE  B1  Contribution  to  the  scintillation  index  (bend-passed  from  0.1 
to  1.5  Hz)  as  a function  of  distance  along  the  line  of  sight 
for  three  source  angular  diameters  (6).  The  calculation  was 
done  for  Fq  « r , q - 3.6,  solar  wind  velocity  V - 400  km/ sec, 
and  34.3  Mix.  The  locations  of  the  centers  of  slabs  7,  13,  19 
and  20  are  indicated. 

FIGURE  B2  (a)  Numerically  calculated  ranges  of  v^.  Each  range  is 

labelled  according  to  source  angular  diameter,  and  they  are 
grouped  according  to  the  choices  of  slabs  to  describe  the 
medium.  Variations  in  the  medium  profile  are  Introduced 
within  each  choice  of  slabs  by  varying  the  elongation, 
resulting  in  a range  of  V3  rather  than  single  points. 

Elongation  was  varied  between  60*  and  120*.  The  solar 
wind  velocity  was  set  at  V ■ 400  km/ sec,  power  lav  slope 
q ■ 2.8. 

(b)  Same  as  (a),  for  q - 3.6. 

FIGURE  B3  Data  from  June-July  1976  for  3C48  and  3C273  presented  in  the 
(a,  v)  plane.  Observations  at  34  MHz  from  the  University  of 
Iowa  G0C0A-Cros8  radio  telescope  and  at  38  MHz  from  the 
University  of  Maryland  TPT  radio  telescope. 


FIGURE  B4  Subset  of  the  data  in  Figure  B3,  adjusted  to  a velocoty  V - 
400  km/sec.  Ranges  of  from  Figure  B2  are  selectively 
Included  for  comparison  with  the  data.  The  3C48  points  are 
consistent  with  a 0. 5"~3.0"  source  with  an  absence  of  tur- 
bulence at  latitudes  > 10°  (hence  large  distances,  >0.7  AU 
along  the  line  of  sight),  and  a 3.0"-7.0"  source  for  3C273 
with  a spherically  symnetric  medium.  The  best  matches  are 
for  q ■ 2.8.  The  unadjusted  v^'s  for  each  source  are  Included 
for  comparison.  Note  that  before  correction  for  velocity,  the 
v^' 8 for  3C48  are  much  more  separated  from  those  for  3C273, 
and  that  only  three  of  the  3C48  points  lie  in  the  range  of  a 
spherically  symmetric  medium  with  V « 400  km  s * while  most 
of  the  3C273  points  lie  in  that  range. 
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These  rainfall  records  were  all  measured  at  the  Anza-Borrego 
Desert  State  Park  headquarters  weather  station  located  approximately 
20  km  from  the  Clark  Lake  Radio  Observatory.  Each  month  in  which  there 
was  more  than  l.S  inches  of  rain  usually  resulted  in  submersion  of  the 
dry  lake  bed  under  a few  inches  of  water.  Monthly  totals  of  more  than 
3 Inches  resulted  in  flooding  to  a depth  of  more  than  6 inches,  and  the 
August,  1977  rainfall  resulted  in  flooding  to  a depth  of  more  than  15 
inches  (Figure  Cl).  Submersion  of  the  d.c.  control  line  cables  for 
extended  periods  of  time  in  the  highly  alkaline  standing  water  resulted 
in  rapid  corrosion  of  the  cables.  Numerous  short-circuits  also  developed, 
many  occurring  several  months  later  as  a result  of  water  retention  and 
continuing  corrosion  inside  the  multiconductor  cable  jacket.  These 
problems  compelled  us  to  install  a new,  open  wire,  elevated  d.c.  control 
line  feed  system. 

We  will  also  note  that  in  March,  1977,  winds  well  in  excess 
of  100  mph  did  extensive  damage  to  the  telescope.  The  subsequent  repair 
and  mechanical  strengthening  of  the  telescope  makes  it  far  more  resistant 
to  such  winds  in  the  future. 
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RAINFALL  TOTALS  BEGINNING  1962 
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